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A COMPUTER PROGRAM FOR FITTING SMOOTH SURFACES 
TO AN AIRCRAFT CONFIGURATION AND OTHER 
THREE-DIMENSIONAL GEOMETRIES 

Charlotte B. Craidon 
Langley Research Center 

SUMMARY 


A digital computer program (D3400) that uses a three-dimensional geometric tech- 
nique for fitting a smooth surface to the component parts of an aircraft configuration is 
presented. The resulting surface equations are useful in performing various kinds of 
calculations in which a three-dimensional mathematical description is necessary. 

Program options may be used to compute information for three -view and ortho- 
graphic projections of the configuration as well as cross-section plots at any orientation 
through the configuration. These operations were implemented to validate the usefulness 
and versatility of the surface equations. Output from this program has been used to drive 
Calcomp, Gerber, and Varian plotters and for on-line display on a cathode-ray-tube device. 

The aircraft (Harris) geometry input section of the program may be easily replaced 
with a surface point description in a different form. Therefore, the program could be of 
use for any three-dimensional surface equations. 

At the present time, the program can only be applied to relatively smooth surfaces; 
that is, there must be no abrupt changes in curvature. This deficiency is overcome to 
some degree by using the airplane component parts or, stated another way, by using a 
collection of surfaces. 


INTRODUCTION 


Aerospace vehicles, automobiles, and ships are examples of objects which require 
smooth curved surfaces to establish their exterior shapes. An integral part of the design, 
development, and manufacture of these objects is the construction of surface models which 
can be analyzed for their interaction with the environment in which they are to operate. 
The most useful models from the point of view of versatility and exactness of definition 
are mathematical models. 

The simplest mathematical model of a three-dimensional surface is a set of planes 
which are defined by points and approximate the curved surface. In order to obtain an 



accurate definition using a discrete set of planes, a large number of points on the surface 
must be defined. Preparing and manipulating the data which yield a planar approximation 
of a surface is laborious if an accurate definition is desired. Another difficulty with 
planar approximation occurs when cross sections or contours of curved surfaces are nec- 
essary. Planar approximation yields a very rough cross section or contour unless an 
extremely large number of points are used to define the surface. 

In recent years a high-order accurate method for mathematical modeling of smooth 
three-dimensional surfaces has been developed. (See ref. 1.) This method is based on 
approximating an arbitrary surface by piecing together surface "patches." Each patch is 
defined by four boundary curves and is bicubic with respect to two parametric variables 
in the interior. A patch is therefore defined by four corner points, the first derivative of 
the corner points with respect to two parametric variables, and the cross derivatives of 
the corner points with respect to the parametric variables. The patch- equation definition 
yields a smooth representation of an arbitrary surface with relatively few points of defini- 
tion. It also yields smooth approximations to cross sections and contour plots. 

The purpose of this report is to describe a computer program which is based on the 
use of sets of bicubic patches to define a relatively smooth surface. (See ref. 2.) In 
particular the data description for the program is oriented toward aircraft configurations. 
This allows the organization of data for the various components to be identical with the 
data used for several standardized aerodynamic analysis computer programs. (See ref. 3.) 
The aircraft data description has become known as the Harris Wave Drag geometry. 

The program can also be used to model arbitrary three-dimensional objects by 
using an alternate data input format. The data -point input to the program is not required 
to be equally spaced in any coordinate variable. However, there are some restrictions on 
the number of points in the descriptive lines for the same surface. 

A three-dimensional parametric cubic spline technique is used to curve fit the input 
data points roughly describing the surface. From the curve fit, the derivatives of the sur- 
face patches with respect to the parametric variable at the corner points are established. 
The cross derivatives of the patch representations with respect to the parametric vari- 
ables are not used in this program. The values of the corner points and the derivatives 
at the corner points constitute the information necessary to solve the patch equation. In 
this way 36 pieces of information are required to define a patch; however, only 12 pieces 
of information must be supplied as input. The remainder is determined from the spline 
fit. Appendix A describes the cubic spline fit technique and appendix B, the patch 
equations. 

The entire aircraft geometry or other three-dimensional object is converted into 
surface patch form. Each patch definition is identical in matrix structure which simpli- 
fies the organization of the computer program and data base. (See appendix C.) All 
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computations, such as rotations, are performed directly on patch equations rather than on 
interpolated x-, y-, and z-coordinates. 

The computer program has the ability to display the orthographic projections of the 
input description of the surface and the enriched description of the surface based on the 
patch definition. (See appendix D.) The desired angles of orientation for viewing the 
surface are inputs to the computer program and the transformation based on these angles 
is applied to the patch definition. An option of the program tests the derivatives normal 
to the surface and deletes those points from the orthographic projection which are facing 
away from the observer. This option gives the program a partial hidden-line capability 
which works very well for convex closed surfaces. Figures 1 to 3 are examples of ortho- 
graphic projections. 

The program is also capable of producing plots of the surface coordinates of a cross 
section at any desired orientation. (See appendix E.) Figures 4 to 6 are examples of 
cross-section plots. The cross-section calculations consist of the simultaneous solution 
of the patch equations and the equation of a plane. The plane is defined by three points 
which are input to the program. 


SYMBOLS 

A,B,C,D parametric cubic spline coefficients 
a,b,c,d plane equation coefficients 

B boundary matrix 

L chord length 

M blending function matrix 

M Mach number 

N surface normal vector 

P a vector whose components are functions of t 

Pj,P 2 ,Pg points used to define a plane 

S component of surface patch equation 
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T diagonal vectors 

t independent variable in cubic equation 

u,w independent variables in patch equation 

V a vector whose components are functions of u and w 

v unit vector 

x,y,z coordinates of a point 

6 pitch angle 

0 roll angle for Mach plane orientation 

0 roll angle 

yaw angle 

PROBLEM DESCRIPTION AND METHOD OF SOLUTION 


The numerical model of the input airplane configuration is assumed to be symmet- 
rical about the XZ-plane (positive Y-side) and may include any combination of components: 
wing, body, pods, fins, and canards. The wing is made up of airfoil sections, the fuselage 
is defined by either circular or arbitrary sections, the pods are defined similar to the 
fuselage, and fins and canards are defined similar to the wings. 

The configuration is usually positioned with its nose at the origin and with the length 
of the body stretching in the positive X-direction. 

The coordinate system used for this program is a right-handed Cartesian coordinate 
system as illustrated in sketch (a). 

Since the modeling technique expects to approximate a smooth surface, sufficient 
input data points with no abrupt changes in curvature should be supplied. A three- 
dimensional parametric cubic spline technique is used for the patch boundary -curve defi- 
nitions in which the coordinates are expressed as cubic functions of one variable. A 
series of adjacent polynomial segments between each given point is used to represent the 
curve. The length of each segment is used as the parameter and later normalized to 1. 
Linear segments are used when a line consists of less than three points. 
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The cubic spline curve fitting subroutine uses a technique from a paper by Timothy 
Johnson of Massachusetts Institute of Technology and is summarized in appendix A. 

The x-, y-, and z-coordinates of a surface patch are each single-valued cubic func- 
tions of two parameters, u and w. The coefficients of these cubics are expressed in 
terms of end points and partial and cross derivatives with respect to the u and w 
parameters. The result is a parametric bicubic representation of three-dimensional 
surfaces. Sketch (b) shows a typical patch. 



Each patch equation requires 48 pieces of information. 

A summary of the bicubic surface patch equation form is given in appendix B and 
the storage file description of the surface patches is given in appendix C. 

The orthographic projections illustrated in this report are created by applying the 
three-dimensional rotation equations directly to the patch equations describing the body 
surface for plotting the aircraft at any desired viewing angle. The rotated patch equa- 
tions are projected into the two-dimensional patch form of the paper plane. An enriched 
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surface may be obtained from the rotated and projected patch equations by holding u 
constant and varying w from 0.0 to 1.0, then by holding w constant and varying u 
from 0.0 to 1.0. 

The orthographic plotting routine also includes a hidden-line option where the nor- 
mal vectors are computed from the rotated and projected form of the patches. A positive 
normal vector indicates that the point is visible and a negative normal vector indicates 
that the vector points away from the viewer and thus is not visible. The method used for 
the orthographic projections is given in appendix D. 

Another routine has been written to compute and plot the surface coordinates of a 
cross section through the body at any desired orientation. The calculations consist of 
the simultaneous solution of the patch equations and the equation of a plane. The method 
is described in appendix E. 


PROGRAM DESCRIPTION 
LABELED COMMON 

The following list contains the FORTRAN variables appearing in labeled COMMON. 


COMMON 

label 

FORTRAN 

variable 

Description 

PATPLT 

XMIN 

Minimum x-value for plotting 


XMAX 

Maximum x-value for plotting 


YMIN 

Minimum y -value for plotting 


YMAX 

Maximum y -value for plotting 


ZMIN 

Minimum z -value for plotting 


ZMAX 

Maximum z -value for plotting 


NOBJ 

Total number of objects (or components) which 
could form an aircraft configuration 

THREED 

ABCDE(8) 

Identification 


HORZ 

X-axis of the paper plane 


VERT 

Y-axis of the paper plane 


TEST1 

Hidden- line option flag 


6 



COMMON 

label 


XSECT 


FORTRAN 

variable 

PHI 

THETA 

PSI 

PLOTS Z 

TYPE 

NOU 

NOW 

ISIDE 

KODE 

ABCDE(8) 

PPL1(3) 

PPL2(3) 


PPL3(3) 

PLOTSZ 

HPAGE 

VPAGE 

INP 

NOU 

NOW 

ISIDE 


Description 


Roll angle 
Pitch angle 
Yaw angle 
Plot frame size 
Type of plot desired 

Number of internal points for each patch in 
u-direction 

Number of internal points for each patch in 
w-direction 

Flag for plotting object or object and its mirror 
image 

Flag for plot-option branch 


Identification 

Origin of cross-section plot and one point in 
three -point-plane definition 

Second point in three-point-plane definition; 
or X-intercept, roll angle, and Mach number 
in plane-angle definition 

Third point in three-point-plane definition 

Scale factor 

Horizontal paper origin 

Twice the vertical paper origin 

Specifies kind of plane input 

Number of points to interpolate for each patch 
in u-direction 

Number of points to interpolate for each patch 
in w-direction 

Flag for examining object or object and its 
mirror image 
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COMMON 

label 


FORTRAN 

variable 


Description 


IPRIN 

Print code 

KODE 

Flag for plot-option branch 

XSTAT 

X-intercept for Mach plane 

THETR 

Roll angle for Mach plane orientation 

XMACH 

Mach number 


OVERLAY ARRANGEMENT 

Program D3400 is set up in the overlay mode and sketch (c) illustrates the overlay 
arrangement. 


Overlay (0,0) 



Sketch (c) 


The control program (0,0) calls in the other parts of the program as they are needed. 
The initialization overlay (1,0) reads cards defining the body surface, converts the input to 
actual units, and temporarily stores the surface description as a series of lines. Cubic 
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spline fairing is performed on curved lines defining the body surface in overlay (2,0) and 
surface patch equations are constructed and temporarily stored. Overlay (3,0) generates 
orthographic plot information using the patch equations, and overlay (4,0) generates cross 
sectional plot information using the patch equations. 
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PROGRAMS AND SUBPROGRAMS 


Program D3400 

Program D3400 (overlay (0,0)) is the control program. This program initiates 
loading and execution of other parts of the program as required. The flow chart and the 
FORTRAN statements for this overlay are as follows: 
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OVERLAY (CSC, 0,0) 

PROGRAM 03 400 ( INPUT=20l ,QUTPUT=201 ,TAPE10=20l, 

1TAP£5= INPUT, TAPE6=0UTPUT , TAP£7i 
C 

C 03400 (SPADE) - SURFACE PATCH OEFINITION EQUATIONS 

C (CONVERTS A SURFACE POINT DESCRIPTION 

C TO THREE DIMENSIONAL SURFACE PATCH EQUATIONS) 

C PRUGRAMtR - CHARCOT T E B. CRAIOON 

C 

COMMON/PAT PL T/ 

IXMIN,XMAX»YMIN»Y MAX »2MIN*ZMAX ,NOB J 
C 
C 

CdC=3LC6C 
RECALL=6HREC ALL 
CALL PSEUDO 
RRi T E (6,20) 

20 FORMAT i 1H I IDX , 2 4H PROGRAM 03400 (SPAOE) - ,34HSURFACE PATCH DEFINI 
I T I ON EQUATIONS///) 

C 

C INPUT SURFACE POINT DESCRIPTION AND PROCESS 

C FUR TAPE 10 AND FUR LABELED COMMON PAT' PL I 

C 

30 CONTINUE 

CALL OVERLAY (CGC, 1,0,0) 

C 

C COMPUTE AND STORE PATCH EQUATIONS 

C 

CALL OVERLAY (COL, 2, 0,01 
C 

40 READ (3,30) I TYPE 

30 FORMAT (14) 

IF (ENUFILE 3 ) 90,60 
oO uO TO (70, BO, JO), I TYPE 

C 

C THREE DIMENSIONAL PLOTS 

C 

70 CONTINUE 

CALL OVERLAY (C3C, 3,0,0) 

GO TO 40 
C 

L CROSS SECTION PLOTS 

C 

BO CALL OVERLAY (C6C,4,0,0) 

GO TO 40 

90 CALL NFRAME 3 CALL CALPLT ( 0. , 0. , 999) $ STOP 

C 

C END OF D 3400 

C 

END 
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Program START 


Program START (overlay (1,0)) reads the configuration description cards and prints 
them, changes the input values to actual units where necessary, computes the minimum 
and maximum dimensions of the given configuration, and temporarily stores the airplane 
description as a series of lines. The flow chart and the FORTRAN statements for this 
overlay are as follows: 
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Read 

XFUS cards 
and 
print 


Write 

gri 

uni 

fin 
d on 
t 10 

— Z 

L 

Search for 
min. and 
max. 

dimensions 

Z 




Change fin 
input to 
actual unit! 

z 

£ 


/Read CANORG, > 
/ XCAN and 
IcANORD cards 
V and print j 


'Read FINORG> 
XFIN , and 
FINORD 
cards 

vand print > 



write pod 
grid on 
unit 10 


Search for 
rain. 

and max. 
dimensions 



/Read PODORG> 
/ XPOD and 
IpODORD cards 
Y and print 


rite body 
grid on 
unit 10 


13 



















o o o 


OVERLAY ( C8C, 1 » 01 
PROGRAM START 

C 

C INPUTS AIRCRAFT SURFACE DESCRIPTION* 

C FORMS INTO DESCRIPTIVE LINES WRITTEN ON TAPE 10, 

C AND COMPUTES MINIMUMS AND MAX1MUMS 

C 

COMMON/PATPLT/ 

1XMIN,XMAX, YMIN. YHAX , ZM I N , ZMAX ,NOB J 
C 

DIMENSION BLOCK (73001 
C 

DIMENSION XAF 1 30 1 , WAFORG ( 20 » 4 ) » W AFORD (20,3 *30) , TZORD ( 20 , 30 1 
tOUI VALENCE ( bLUCK * XAF i , ( BLOCK! 31 1 , WAFORG)-, 

1 ( BLOCK! Ill), wAFU.KD) , (6L0CK(19il) , TZORD) 

C 

DIMENSION XFUS(30,4) ,ZFUS(30,4) ,FUSARD( 30,4) ,FOSRAD( 30 ,4) , 

1 SFUS ( 30,30 , B ) 

EQUIVALENCE (BLJCK.XFUSJ , (BLOCK( 121) ,ZFUS) , (BLOCK (241) .FUSARD) , 

1 ( BLOCK! 361 ) , FUSR AD) , ( BLUCK ( 241 ) ,SFOS) 

C 

DIMENSION PODURG (9 , 3 ) , XP OD (9 , 30 ) ,P0D0RD(9, 30) , XP0D1 (9, 30) 
EQUIVALENCE (BLOCK , PODURG) , ( BLOCK! 2b) »XPOD) , (BL0CK(298) .PUOORO), 

1 ( BLOCK (5o8) , XP0U1) 

C 

DIMENSION FI NORG (6,2,4) ,XFIN(6,10) , FIN0RD1.6, 2, 10) , 

1FINX2(6,2, 10 ) , F INX3( 6,2, 10) 

EQUIVALENCE l BLOCK, F I NORG) , (BLOCK (49) ,XFIN), (BLOCK! 109) ,FINORD) , 

1 ( BLOCK (229),FINX2) , ( BLUCK ( 349) ,F INX3) 

C 

DIMENSION CANORG (2,2,4) ,XCAN(2,I0) ,C ANORD (2 , 2, 10) , 

1 CAnORI ( 2 , 2 , 1 0) , C ANURX ( 2,2,10) 

equivalence ( block, c anorg ) , i bluck( it ) ,xcaw) , (block(37) ,canord) , 

1 (BLOCK! 77) , C ANOR 1) ,( BLuC K ( 1 1 7 ), C ANURX ) 

C 

DIMENSION ABC l B ) , ABC 0 ( 8 ) ,ANSIN(30) , AN COS ( 3 0) » NAME ( 2 ) 

DIMENSION NRADX(9) ,NF0RX(4) 

DIMENSION ALRT (31,3) 

DATA PI/3. 19159265/ 

C 

REWIND 10 

10 FORMAT ( 8A1 0 ) 

20 FORMAT (1XBA10) 

30 FURMAT (10F7.0) 

READ ID CARD AND CARD OF CONTROL INTEGERS 

READ (5,10) ABC 
IF (ENDFILE 5) 35,40 
35 CALL NFRAME $ CALL C ALP L T ( 0. ,0. , 999 ) $ STOP 
40 CONTINUE 

WRITE (6,50) ABC 

50 FURMAT (23X, 34HAIRCRAFT CONFIGURATION DES0R1 PTIUN//1X8A10/) 

READ (5,10) ABCD 
WRITE (6,60) ABCD 
oO FURMAT (1X8A10/) 

DECODE ( 72, 7 0, ABCD) JO, J 1 , J 2 , J 3, J4, J5 • J6,NWAF , NWAFOR , NFUS , ( NRADX ( I 
l ) ,NFORX( I ) ,1=1,4) » NP » NPODQR , NF» NF l NOR, N CAN ,NCANOR 
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70 FORMAT (2413) 

C 

WRITE (10) ABC 
N08J=0 
C 

C REFERENC E AREA 

C 

IF (JO.NE.l) GQ TO 80 
READ (5,10) AGCO 
WRITE (6,20) ABCO 
C 

C WING 

C 

80 J J=1 A8S( J l ) 

IF (JJ.NE.l) GO TO 290 
N=IABS(NwAFOR) 

NREC=(N+9) /10 

1 l=-9 

12=0 

DO 90 NN=l ,NREC 
READ (5,10) ABCD 
WRITE (6,20) ABCO 
11=1 1+10 
12=12+10 

DECODE (70,30, ABCD) ( XAF ( I ) , I =1 1 , 12 ) 

90 CONTINUE 

DO 100 1 =1 »NWAF 
READ (5,10) A6CD 
WRITE (6,20) ABCO 

DECODE (26 , 30, ABCD) ( WAFCRGI I , J ) , J= l , 4) 
100 CONTINUE 

IF (Jl.LT.O) GU TO 130 
DO 120 NN= 1 , NWAF 
1 1 =-9 
12=0 

DO 110 N 1= 1 , NREC 
READ (5,10) AoCD 
WRITE (6,20) AbcD 
11=1 1+10 
12=12+10 

DECODE (70, 30, ABCD) (TAORD(NN»I) ,1=11,12) 
110 CONTINUE 
120 CONTINUE 
GO TO 150 

130 UU 1 AO 1 = 1, NWAF 
DO 140 K= 1 , N 
140 TZURD1 I ,K j=0. 

150 L=1 

IF (NwAFOk.CT.O) L=2 
DO 170 NN=l»NwAF 
DO 170 K=i,L 
1 1=— 9 
( 2=0 

DU 160 Nl= 1 » NREC 
READ (5,10) ABCD 
WRITE (6,20) ABCD 
1 1=11+10 
12=1 2+10 
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160 

170 


iao 

1*0 


c 

c 

c 


200 

210 


220 

2JU 

C 

c 

c 


c 

c 

c 


OECGUE (70fiQ#AbCUJ (WAFGKlHNN 

CUNT I NUE 

CONTINUE 

IF { NWAFUR.LT. Oj GO TO 190 
Ou 180 NN- 1 y NttAF 
00 180 K=i t N 

rtAFORG(NNf 2f K AFGRO( N N » 1 , K ) 

continue 

NWAF0R=IA8S<NwAFUR) 

N*=N*AFGR 
Jl=l A8S( J1 J 




CHANGE TO ACTUAL UNITS, COMPUTE MIN1MUMS 


00 210 1=1 »NWAF 
E= . 0 1#NAFURG ( I ,4) 
E3=W AFURG ( 1,3) 


Ou 200 0=1 ,nwafur 

WAFOROI I , 1 ,J )=E*WAFORD( I , 1 , J J +£ 3+TZORDI I 
WAFOROI I,2,J)=— E*WAF0RU(I,2,J 1+E3+TZ0R0I 
WAF0R0(I,3,J) = »aF0RG(I,U + E*XAF{J) 

CONTINUE 


.J) 
1, J) 


XMIN=XMAX=w AFORG 11, 1 J 
YMAX=WAFORG< 1,2) 

YMIN=WAFURG( 1,2) 

ZM IN = ZMAX= wAFQRO ( 1 ,1,1) 

00 230 N=1 »NWAF 

XMAX=AHAXi UMAX, wAFQROl N,3,NW)) 
XMIN=AMIN1 (XMIN, WAFOROI N, 3, 1)1 
YMAX= AMAXi ( YMAX, WAFuRGl N»2) ) 
YMIN=AMIN1 ( Y Ml N, WAFORGI N, 2) ) 

DO 220 NN= 1 , Nto 

ZMAX= AMAX1 I ZMAX , WAFORQ ( N , i , NN) ) 

ZMIN-AMINKZMIN, WAFOROI N » 2 »NNJ ) 

CUNTINUE 

CONT INUE 


AND MAXIMUMS 


WRITc LINE TAPE 


NN=2 


NCOMP= l$NAME I 1 ) = 10HW ING *NAMEI2)=1 

WRITE (10) NN, NGOMP, NAME *NN , 

N0BJ=NUBJ«-1 


00 280 1=1,2 

WRITE (10) NWAF,NWAFGR,NN,NN,NN 
RKK=( I-i)*lNWAFOR*i) 

KK=(-1 )**( i + ii 


SETUP SPANWISE LINES 


240 

250 


OU 250 K=1 ,NWAFUR 
NN=KKK+KX*K 


00 240 N=1,NWAF 

ALRT IN»1)=WAF0R0(N,3»NN) 

ALRT I N,2)=wAF0RG(N,2 ) 

ALRT IN»3)=WAF0RD(N» I ,NN) 

CONTINUE 

WRITE (10) ( ( ALRT (N»N3 ) , N= l, NWAF ) 
CONTINUE 


»N3=l,3) 
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SETUP ST RE AMWi St LINES 


C 
C 

c 

DU 2 TO NN=l»NwAF 
DU 260 K=1 , NwAFUk 
N=KRR+KK*X 

ALRT (K»i)=WAF0RDlNN*3»N) 

ALKT (K,2)=WAFoRG(NN,2) 

ALRT (K, 3) =WA FORD INN, l,N) 

260 CUNT INUE 

WRITE (10) I (ALRT(N,N3) ,N=l,NW) ,N3=1, 3) 
270 CONTINUE 
280 CONTINUE 
C 

C FUSELAGE 

C 

290 JJ=I ABSt J2 ) 

IF l JJ.NE. I) GU TO 590 
J2TEST=3 

IF (U2.EQ.-1 .AND.Jo.EU.-U J2IEST=1 
IF ( J2.EU.-1.ANU.J6.E0.0) J2TEST=2 
IF l J6.tw. I) J2TEST= l 
J2=l 

DO 410 NFU— I »NFUS 
NRA0=NRADX(NFU) 

NFUSGR=NFORX (NFU ) 

N=NFUSUR 
NREC=(N+9)/10 
1 l=-9 
12=0 

OU 300 Nl= I » NREC 
READ I 5 » LO ) ABCD 
WRITE 16,20) AbCD 
11 = 11+10 
12 = 12+10 

DECODE (70 ,30, A3CD) I XFUS I I , NFU ) , 1= 1 1 » I 2) 
300 CONTINUE 

IF ( J2TEST .NE.2) GO TO 320 
1 1=-9 
1 2=0 

DO 310 Nl= l , NREC 
READ (5,10) ABCD 
WRITE (6,20) A6CD 
11 = 11+10 
12 = 12+10 

DECODE (70,30, ABCD ) ( 2FUS ( It NFU ) ,1=11,12) 
310 CONTINUE 
GU TO 340 

320 DO 330 1 = 1, N 
330 2FUSI l,NFU)=0. 

340 IF 1 J2TEST.NE.3) GO TO 380 
NCARO= (NRA0+9) /I 0 
DO 370 LN=I» N 
DO 360 K=1 ,2 
KK=K+(NFU-1) *2 
11 = 10 
1 1=— 9 
12=0 
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DO 350 NN= l , NCARD 

IF (NN.EO. NCARD) 1 I = MOU (NR AD, 10) 

IF ( I L.EQ.O) 11=10 
11 = 11+10 
12 = 12+11 

REAO (5,10) ABCO 
(•RITE (6,20) A BCD 

DECODE (70, 30, AB CD) ( SFUS( I , LN,KK) , 1=11,12) 
350 CONTINUE 
360 CONTINUE 
570 CONTINUE 
GO TO 410 
3d0 I 1=— 9 
12=0 

DO 390 Ni= l , NR EC 
READ (5,10) ABCD 
WRITE (6,20) ABCD 
11 = 11+10 
12 = 12+10 

DECODE (70, 30, ABCO) ( FU S ARD ( I, NFU), 1=11,12) 
390 CONTINUE 

DO 400 1=1, N 

400 FUSRADII *NFU)=SwRT(FUSARD( I , NFU) /PI ) 

C 

410 CONTINUE 
C 

C fuselage min and max 

c 

IF ( Jl.Nc.O) GO TO 430 
XMIN=XFUS( 1,1) 

XMAX=XFUS( 1, 1) 

IF ( J2TEST.EQ.3) GO TO 420 
YMIN=FUSRAD(1*1) 

YMAX=FUSRAD( 1,1) 

ZMIN=— FUSRAD (1,1 ) + ZFUS ( 1, l) 

2MAX=FUSRAD( 1,1 ) +ZFUS (1,1) 

GO TO 430 

*,20 YMAX=SFUS( 1,1,1) 

YMIN=SFUS( 1,1,1) 

ZMIN=SFU5( 1, 1,2) 

ZMAX=SFUSI 1, 1,2) 

+30 DU 470 N=i ,NFUS 
NkAD=NRADX(N) 

NFUSOR=NFORX(N) 

aMIN=AMIN1(XMIN,XFUS(1,N)) 

XMAa=AMAX1 ( XMAX, XFUS(NFUSGR,N) ) 

DO 4o0 NN= 1 , NFUSOR 

IF ( J2TEST.E0.3) GO TO 440 

YMAX=AMAX1 ( YMaX, FUSRAU ( NN ,N ) ) 

YM IN = AM IN 1 ( YMIN , FU3R AD ( NN »N ) ) 

2MAX=AMAX1 (2 MAX, FUSRADt NN»N) +ZFUS(NN,N) ). 

2M IN=AMI N1 ( EMIN, —FU SR AO ( NN* N ) +ZFUS ( NN ,N ) ) 

GO TO 460 
440 KK= 1 +1 N-l ) *2 

DO 450 NR= 1 , NR AO 

YMIN=AMIN1 (YMIN,SFUS(NR,NN,KK)) 

YMAX=AMAX1 ( YMAX, SFUS ( NR , NN »KK) ) 
ZMIN=AM1N1(ZM1n, SFUS (NR, NN,KK+ i) ) 

430 ZMAX = AMAX1(ZMAX» SFUS (NR,NN,KK+l) ) 
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460 CONTINUE 
470 CONTINUE 
C 

C WRITE LINE TAPE 

C 

JJN=2*Nl=l$NAME(l)=lOHFUSELAGE $NAME ( 2 1=1. OH 

NOB J — NOB J + NF US 

DO 580 NFU=1 »NFUS 

NRAO— NKADX INFUi 

NFUSOR=NFORX (NFU ) 

WRITE (10) N i t J JNt NAME > N 1 1 NI 
WRITE (10) NFUSUR, NRAU,N1,N1,N1 
NAN=NRAD 

IF ( J2TEST . t w. 3 J GO TO 400 

FANG=(NRA0-1)*2 

DELE = 6. 2 63 1 8 53 /F ANG 

DO 480 N=I * NAN 

£=N~ 1 

ANS IN( N)=S IN l E*UEL£*4. 7 12389) 

480 ANCOS ( N)= CQSIE^U EL E+4. 712389) 

490 CONTINUE 

KK=l*(NFU-l)*2 

C 

C SETUP STREAMWISE LINES 

C 

DO 530 N=1 t NAN 
DO 520 NN= 1 • NFUSUR 
ALRT (NN» 1)=XFUS( NN»NFU) 

IF ( J2TEST.EU.3) GO TO 500 

ALRT ( NN ,2 ) =F USRAD( NN,NFU)*ANCUS(N) 

ALRT (NN,3)=FUSRA0(NN ,NFU ) *ANS IN { N ) +ZF US INN ,NFU) 
GO TO 510 

500 ALRT(NN,2)=SFUS(N,NN,KM 

ALRT (NN»3)=SFUS(N,NN,KK+1) 

510 CONTINUE 
520 CONTINUE 

WRITE (10) I (ALRT(N,N3) , N= 1, NFUSOR) ,N3= 1, 3) 

530 CONTINUE 
C 

C SETUP LINES AROUND BODY 

C 

DU 570 N=i, NFUSUR 

DU 560 NN= l « NAN 

ALRT (NN, 1)=XFUS(N,NFU) 

IF ( J2 TEST . EQ.3) GO TO 540 

ALRT ( NN >2 ) =FUSRAo( N» NFU ) *ANCOS( NN) 

ALRT(NN,3)=FUSkA0IN,NFU)*ANSIN(NN)+ZFUSIN,NFU) 

GO TO 550 

540 ALRT(NN,2)=SFUS(NN,N,KK) 

ALRT ( NNt 3) =SFUS I NN »N » KK + 1 ) 

550 CONTINUE 
560 CONTINUE 

WRITE (10) ( I ALRT ( N, N3 ) , N=1 » NAN) » N3=l ,3 ) 

570 CONTINUE 
580 CONTINUE 
C 

C NACELLES 

C 
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5*0 CONTINUE 

IF (Ji.NE.l) 00 TO 730 

N=NPGOOR 

NREC=(N+*)/10 

OU 620 NN= 1 » NP 

READ (5,101 AbCD 

WRITE (6,20) AbCU 

UECOGE (21,30, A3CD) ( PUOORG ( NN, 11,1-1,3) 

1 1=- 9 
12-0 

DJ 600 Ni=l,NREC 
READ (5,101 ABCD 
WRITE (6,20) AbCU 
1 1 = 11*10 
12 = 12+10 

DECODE ( 70 , 30, ABCU 1 ( XPOO ( NN , 1 1 , 1= 1 1 , 1 2 1 
600 CONTINUE 
Il=-0 
12=0 

DO 610 Nl= 1 » NREC 
READ (5,10) A3C0 
WRITE (6,20) AbCD 
( 1 = 11+10 
12 = 12+10 

DECODE ( 70 ,30, ABCD) 1 PODCRD ( NN, i ) , 1=1 1 , (2) 
610 CONTINUE 
620 CONT INUE 
C 

C COMPUTE ACTUAL X , MI N I MUM ,MAX 1MUM 

C 

DO 630 N=l,NP 
DO 630 NN= 1 , NPODOR 

630 XPOD1(N,NN)=XPOD(N,NN)+PUDORG(N, 1) 

IE ( JUNE. O.OR.J2.NE.O) GO TO b40 
XMI N=XP0D1 (1,1) 

XMAX=XP001 (1 , NPODOR) 
YMIN=P00URG(i,2)+P0G0RD( 1,1) 

YMAX= PODORGI 1,2) +PUD0R0 (1,1) 

ZM1N=PODORG ( 1,3)— PODORO( 1,1) 

ZMAX=PODORG( 1 , 3 ) +POOORO ( 1 , 1 ) 

640 DO 660 N= 1 , NP 

XM1N= AMINi (XMIN, XPOO 1 ( N , 1)) 

XMAX=AMAX1 (XMAX, XP0D1 (N, NPODOR) ) 

DO 650 NN=1, NPODOR 

YMIN= AMINI (YMIN, PUOORU ( N , NN ) +P000RG( N, 2 ) ) 
YMAX=AMAXl(YMAX, PODORDI N ,NN ) + PODORGI N ,2 1) 
ZMIN=AMI N1 (ZMIN, P00URG(N,3 )— PQDORD ( N, NN ) ) 
650 ZMAX= AMAX1 ( ZMAX , PODORGI N , 3 l+PODURD ( N , NN ) ) 

660 CONTINUE 

DATA NAN2/4/ ,P IPL/4. 712 389/ 

NANG l=NAN2+l 

NANG2=2*NAN2+1 

F ANG=NAN2*2 

DELE=o .2831 853/F ANG 

DO 67 0 N=i , N ANG2 

E=N— 1 

EE=E*DELE 

ANSI N( N)=SIN(tE+PIPL) 

670 ANCOS ( N) =COS (EE+P1PL) 
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WRITE LINE TAPE 


L 
C 
C 

J JN=3*NAME ( 1 ) = lOHPQDS 4>NAME( 2 ) = 10H 

NOaj=NObJ+NP 
DU 720 NPl=l ,NP 
1-2 

IP (POUuRG(NPl,2).EQ.O.) 1=1 
WRITE (10) l * JJN.NAME, I , I 
DU 7 20 1=1,2 

IP ( l.Eu.2.AND.P00UkG(NPl,2) .Eg.O.) GO TO 720 
WRITE (10) NPOOOR, NANG l ,1,1,1 
C 

C SETUP STREAMWISE LINES 

C 

DO 690 K=1,NANG1 
NN= I I— l) *NAN2*K 
DU 600 N=1 » NPOOOR 

ALRT (N, l)=XPOO(NPl »N) +PUDORG( NP 1,1) 

ALRT ( N ,2 )=POOURO ( NPl , N) *ANCOS(NN) *PGDORG(NPl ,2) 

ALRT (N,3) = P00UK0 (NPl , N) * ANS I N(NN ) NPODORGI NPl , 3 ) 
o80 CONTINUE 

WRITE (10) ((ALRT(N,N3) ,N=1, NPOOOR), N3=l, 3) 

690 CONTINUE 
C 

C SETUP LINES AROUND PODS 

C 

DU 710 N=1 tNPODUR 
DU 700 K=1 ,NAnG1 
NN=( I — 1 ) *N AN2+K. 

ALRT ( R » 1 ) = XPOOI NPl »N) +PG0GRG (NP 1 ,1) 

ALRT (K,2)=PUDGR0 (NPl ,N) *ANCGS (NN ) +P0D0RG( NP 1 , 2 ) 
ALRT(K,3)=P0D0RD(NP1 ,N) * ANS I N ( NN ) +P0D0RG1NP1 , 3 ) 

700 CONTINUE 

WRITE (10) ( ( ALRT ( K , N3 ) , K=1 , NANG 1) »N3= l » 3) 

710 CONTINUE 
720 CONTINUE 
C 

C PINS 

C 

730 CONTINUE 

IF ( J9.NE. 1) GO TO 890 
N=NF I NOR 
DO 790 NN= 1 , NF 
HEAD (5,10) A8CU 
WRITE (6,20) AbCD 

DECODE (56, 30, AbCD) ( ( F I NORG ( NN , I , J ) , J= 1 , 9 ) , 1= 1, 2 ) 

READ (5,10) AbCD 
WRITE (6,20) AbCD 

DECODE (TO, 30, AbCD) ( XF IN ( NN , I ) , 1= 1 , N) 

READ (5,10) AbCD 
WRITE (6,20) AtiCD 

DECODE (70, 30, AbCD) ( F INORU(NN, 1 , J) , J=1 ,N) 

790 CONTINUE 
C 

C CHANGE TO ACTUAL UNITS, COMPUTE MINIMUMS AND MAXIMUMS 

C 
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00 760 LU=l,NF 
DO 7<>0 1 = 1,2 
J=3— I 

E=.01*F1NURG(LU, J,4) 

E2=F INORGI LQ , J ,2.) 

00 750 K=1 »NFIN0R 
EE=F INORO ( LQ « 1 , K) *£ 

FINQRDILO, J, K)=E2-*-EE 
FINX2(LQ,J,K)=E2-EE 

750 F1NX3(LU,J,K)=FIN0RG(LQ, J, 1)+E*XFIN(LQ,K) 

760 CONTINUE 
C 

IF ( JI.NE.O.OR.J2.NE.O.OR.J3.NE.O) GO TO 770 
XM IN=F I.NOKG ( 1,1,1) 

XMAX=F INORG ( 1,1, 1) 

YM1N=F INORGi 1, 1,2) 

YMAX=F INORG ( 1,1,2) 

ZM 1 N=F 1N0RG (1,1,3) 

ZMAX=F INORGI 1,1,3) 

770 00 780 N=1 ,NF 

ZMIN=AMIN1(ZMIN,FIN0RG(N,1,3) ) 

ZMAX= AMAX1 (ZMAX, FIN0RG1 N ,2,3 J ) 

00 7 60 N2=i , 2 

XMIN=AMIN1 (XMIN,FINCRG(N,N2,1)) 

XMAX= AMAX1 (XNAX, FI NX 3(N , N2 , NF INOR)) 

00 780 NN=1,NFINUR 
YMIN=AMIN1(YMIN»FINX2(N*N2»NN)) 

YMAX=AMAX1 ( YMAX, FI N0R0( N ,N2,NN) ) 

780 CONTINUE 
C 

C WRITE LINE TAPE 

C 

JJN=4*NAME(1 )=10HFINS $NAME(2)=10H 

NU6J=N06 J+NF 

NK2=2 

00 880 NF1=1,NF 
1=2 

IF (FIN0RG(NFl,l,2).EQ.0.) 1=1 
WRITE (10) I » J JN »NAME ,1,1 
DO 8 70 NN2=i ,2 

IF (NN2.Ey.2.ANU.FlNOkG(NFl,i,2).EQ.O.) GO TO 870 
WHITE (10) NFI NOR, NK2, 1, 1,1 
11=1 
12=2 

IF (NN2.EU.1) GO TO 790 
1 1 = 2 * 12=1 
790 CONTINUE 
C 

C SETUP HORIZONTAL LINES 

C 

00 810 N=l , NF INOR 
ALRT(N,1)=F1NX3 ( nFI, 1 1,N) 

ALRT ( N,3) =F INORO (NFI ,11,3) 

IF (NN2.EU.2) GO TO 800 
ALRT (N,2)=Fi NORO l NFI , Il,N) 

GU TO 810 

800 ALRT (N,2)=FI NX2 (NF1,11,N) 

810 CONTINUE 


22 



WRITE (10) ( (ALRT(N»N3) ,N=l,NFIN0R) ,N3=1,3) 

DO 830 N=1 1 NF1N0R 

ALRT (N,i)=FI NX3( NF 1 » 12, N) 

ALRT ( N » 3) =FI NORG ( NF1 , 12,3) 

IF (NN2.cW.2I GO TO 820 
ALRT ( N , 2 ) =F l NGRD (NF1 , 12 , N ) 

GO TO 830 

820 ALRT (N»2)=FI NX2 (NF1»I2»N) 

830 CONTINUE 

WRITE (10) ( (ALRT ( N , N3 ) » N=1 » NF INQR) »N3=i» 3) 

C 

C SETUP VERTICAL LINES 

C 

DO 860 NN=1, NFINUR 

ALRT (l,l)=FINX3(NFl, U.NNJ 

ALRT (2,l)=FINX3( NF 1 , 12, NN) 

ALRT (1,3) =FI NORG (NF1 , 11,3) 

ALRT ( 2, 3) =FI NORG (NF1 , 12,3) 

IF (NN2.EQ.2) GO TO 840 
ALRTU,2)=FIN0R0(NFl, ll,NN) 

ALRT (2»2)=FINQR0(NFl»I2»NN) 

GO TO 850 

840 ALRT I 1 » 2 ) =F1 NX2( NF 1 , II, NN) 

ALRT (2,2)=FINX2(NF1, 12, NN) 

850 WRITE (10) ( (ALRT(N,N3) ,N=l, 2) »N3=1 ,3) 

860 CONTINUE 
870 CONTINUE 
880 CONTINUE 
C 

C CANARDS 

C 

890 CONTINUE 

IF (J5.NE.1) GU TO 1080 
N=IA8S(NCAN0R) 

DO 920 NN=1,NCAN 
READ (5,10) A8CD 
WRITE (6,20) AbCD 

DECODE (36, 30, AbCD) ( l C ANORG ( NN , I , J ) , J= 1 ,4 ) , 1= 1 , 2 ) 

READ (5,10) A6C0 
WRITE (6,20) AbCD 

DECuDE (70, 30, AbCD) (XCAN(NN, I) , 1=1, N) 

READ (5,10) AbCD 
WRITE (6,20) AbCD 

DECODE (70, 30, AbCD) (CANORD( NN, 1 , J) , 0=1 »N) 

IF ( NCANUR . L T. 0) GU TO 910 
DO 900 J=1,N 

900 C ANUR 1 (NN, 1 , J) =C ANORD( NN , 1 , J ) 

GO TO 920 

910 READ (5,10) ABCD 
WRITE (6,20) ABCD 

uECUOE ( 70, 30, AbCD) (CANORKNN, 1 , J) , J=l ,N) 

920 CONTINUE 

NCANOR= I ABSi NCANUR) 

NC=NCANOR 

C 

C CHANGE TO ACTUAL UNITS, COMPUTE MIN1MUMS AND MAXIMOMS 

C 
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00 950 NN= 1 , NCAiX 

00 940 K=l,2 

1 = 3-K 

E=.01*CAN(JRG(NN, 1,4) 

E3=CANQRG(NN, 1,3) 

DO 930 J = 1 » NCANOR 

CANOROINN, I , J) = £ *C ANORO INN, 1 , J) +E3 
CAN0K1 ( NM * I; J )=-E*CANORl(NN, 1,J).+E3 
930 CANORX(NN, l » J ) =C AMOKG ( NN t I . i» «■£* XCANI NN , J ), 

940 CUNTINUE 
950 CONTINUE 

IE ( JI.NE.O.OR.J2.NE.O.OK.J3.NE.O.OR. J4.N6.0) GO TO 960 
XMIN=CANORX( 1, i, 1) 

XMAX=CANURX( 1 fit NCANOR i 
YMIN=CAN0RG(1,2,2) 

YMAX=C ANURGI 1,2,2) 

ZMIN-CANOR 1 ( 1,1,1) 

ZMAX=C ANORO (1,1, 1) 

960 00 990 NCA=1,NCAN 

YMIN-AMIU1 ( YM1N, CANORGliXCA, l ,2) ) 

YMAX=AMAXl (YMAX, CANORGI NCA,2,2) ) 

00 980 1X2= 1,2 

XM1N=AMIN1 ( XMI N, CANGRX ( NC A » N 2, 1 ) ) 
XMAX=AMAX1(XMAX,CANGRX( NCA,N2, NCANOR) ) 

00 9 70 NN=1, NCANOR 

ZMIN=ANINI(ZMIN,CAN0R1(NCA,N2,NN)) 

970 ZMAX= AMAX1 (ZMAX, CANORDl NCA,N2,NN) ) 

980 CONTINUE 
990 CONTINUE 
C 

C WRITE LINE TAPE 

C 

JJN=5$NAME ( l )= 10HC ANAROS *NAME(2)=10H 

N08J=N0BJ+NCAN 

NK2=2 

DO 1070 NC A= 1 , NC AN 
WRITE (10) NK2,JJN,NAME,NK2,NK2 
00 1060 1=1,2 

WRITE (10) NK2,NC,NR2,NK2,NK2 
KKK=(I-1)*(NC+1) 

KK=(-1)**( H-l) 

C 

C SETUP SPANWISE LINES 

C 

00 1020 K= 1 , NC 
NN=KKK+KK*K 
00 1010 N2 = l , 2 

ALRT (N2, 1)=CANURX( NCA,N2,NN) 

ALRT (N2,2)=CAN0RG(NCA,N2,2) 

IE ( I.Eu.2) GO TO 1000 
ALRT (N2,3)=CANORO( NCA.N2 ,NN) 

GO TO 1010 

1000 ALRT ( N2 , 3 ) =C ANOR 1(NCA,N2,NN) 

1010 CONTINUE 

WRITE (10) ( (ALRT(N2,N3) ,N2=1,2) ,N3=1 ,3) 

1020 CONTINUE 
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c 

c 

c 


SETUP Twu ST REAMWISE LlNcS 


DO 1050 112=1.2 
DO 1040 N=1.NC 
J=KKK+KK*N 

ALRT tN,l)=CAN0RX(NCA,N2, J) 
ALRT(N,2)=CAN0RGINCA,N2,2I 
IF (I.EU.2J GO TO 1030 
ACRT IN,3)=CAN0R0(NCA,N2 » J) 

GO TO 1040 

1030 ALRT ( N . 3 J =C ANUR 1 (NCA »N2 . J J 
1040 CONTINUE 

WRITE ( 1UJ ( ( ALR TIN.N31 » N= l » NCI . N3=i » 3 ) 
1050 CONTINUE 
1060 CONTINUE 
1070 CONTINUE 
1080 CONTINUE 
RETURN 
C 

C END OF START 

C 

END 


! 
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Program SURF 


Program SURF (overlay (2,0)) is the control program for constructing surface patch 
equations. The flow chart and the FORTRAN statements for this overlay are as follows: 


C 


Overlay 

Program 

(2,0) 
i SURF 

S 

7 

/ Write \ 
/identification , \ 
1 min. and max. I 
\ values on 1 
\ unit 7 / 




Write 
control 
integers 
on unit 
l 


call 

PACH 


Compute 
patches 
Kand store , 


) 


C 


D 


OVERLAY (CBC»2fOJ 
PROGRAM SURF 
C 

C CALLS A SUBROUTINE TO COMPUTE PATCHES 

C AND CONTROLS WRITING OF PATCH TAPE 

C 

COMMON/ PAT PLT/ 

1XMIN,XMAX, YM IN, Y MAX t ZM IN, ZMAX.NUdJ 
DIMENSION ABC(dl 
REWIND 7 
REWIND 10 
READ (101 ABC 

CALL RECOUT ( 7 ,2 , 0 ,ABC , 1 » d , 1 J 

CALL RECOUT (7,i,0,XMIN,XMAX,YMiN,YMAX,ZMLN,XMAX,N06Jl 
DO 20 1 = 1 f NOB J 

READ (101 NSURF, J2,J3, J4, J5, J6 
CALL RECOUT ( 7 , 1 , 0 , NSURF , J2 , J3 , J4 , Ji , J6 1 
DU 10 N=i, NSURF 
READ (101 NC0L,NRGW,N3,N4,NS 
CALL PACH INCUL,NR0W,N3,N4,N51 
10 CONTINUE 

20 CONTINUE 

END FILE 7 

RETURN 

END 
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Subroutine PACH 


Subroutine PACH computes surface patch equations from the given grid information 
describing a surface. The patch equations are stored for further use. The flow chart 
and the FORTRAN statements for this subroutine are as follows: 
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SUBROUTINE PACH l NLD , NL S r L 1 , L2, L3 ) 

C 

C CONSTRUCTS SURf-ACc PATCHES WITH THE COMPONENTS 

C EXPRESSED AS CUBIC FUNCTIONS OF TWO PARAMETERS IU AND W) 

C AND WRITES ON TAPE 

C 

DIMENSION PATCH! 4,4,3) , COEF1! 31 , A , 3) , C0EF2! 3 1 , A , 3) 

DIMENSION SLOPE! 31,31,3) , XMAT! A, A) , AL INEl 31 , 3) , ELEN! 3 1 ) , P AT t A , A) 
DATA (XMaT(I)»I=1,1o)/ 2.»-3.,0.» I# ,*2, , 3« ,'0* , 0«i 
il.,-2.»l.,0. *l.»-l.»0.»0./ 

DATA MAXN/31/, EPS/. 00001/ 

NI=NLD— 1 
N2-NLS-1 

CALL RECuUT ! 7 , 1 , 0 , N 1 , N2 , L 1 , L2 , L 3) 

C 

C COMPUTE PARAMETRIC SLOPES IN W DIRECTION 

c 

DO 70 N-l.NLS 

READ (10) ( ( AL I NE ( NN » N3 ) » NN= 1 » NLD ) , N3= 1 ,3 ). 

C 

C CHECK IF CUBIC FAIRING POSSIBLE 

C 

IF ( NLD.LT .3) GO TO 20 
C 

C CHECK FOR A POINT 

C 

Nl=NL0-1 
DO 10 NN=1 ,NL 

El-AB5(ALlNt(NN» l ) -AL I NE ( NN+ 1 , l ) ) 

E2=AbS(ALiNE(NN,2)-ALlNE(NN«-l,2) ) 

£ 3=ABS ( ALI NE ( NN » 3 ) -A L l NE l NN+ 1 ,3 ) ) 

IF ( E1+E2+E3.LE. EPS) GO TO 20 
10 CONTINUE 

GO TO AO 
C 

C COMPUTE LINEAR SLOPES 

C 

20 DO 30 NN=1 , NLD 

DO 30 N3— 1 , 3 

SL0PE!NN,N,N3)=AL1NE12,N3)-ALINE(1,N3) 

30 CONTINUE 

GO TO 70 
C 

C COMPUTE CUBIC SPLINE SLOPES 

C 

AO CALL SPFIT ( MAXN ,NLD, AL I NE, E LEN, COEF 1 , 1 l , 0 , 1 2 , EK , CP , 1 3 , I A) 

NL=NLD-1 
DO 50 NN=l,NL 
DO 50 N3- 1,3 

SLOP E!NN,N,N3)=C0EF1(NN,3,N3) 

50 CONTINUE 

DO 60 N3=l , 3 

SL0PE(NL0,N,N3)=3.*CQcFl INLD-l, l,N3)+2. *COEF 1 ( NL D- 1 , 2 , N3 )*COEP l ( NL 
ID— 1, 3,N3) 

60 CONTINUE 

70 CONTINUE 

C 
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c COMPUTE PARAMETRIC SLOPES IN U DIRECTION 

C FORM PATCHES AND WRITE TAPE 

C 

RcAU l 101 ( ( AL INE( NN tN3 1 »NN= 1 »NLS) t N3- l t3L 

C 

C CHECK IP CUbIC PAIRING POSSIBLE 

C 

IF INLS.LT. 31 GO TO 90 
C 

C CHECK FOR A POINT 

C 

NL=NL3-1 
00 BO NN= I » NL 

EI-ABS ( AL 1 Nc ( NNt 1 1 -AL l NE I NN* 1 » 1 1 I 
t2=ABS( ALINE INN, 2) -A LINE! NN+lt2J) 

E3=A6S(ALiNE(NN t 31-ALlNE( NN+1,3) ) 

IP (E1+E2+E3.LE.EPS1 GO TO 90 
80 CONT INUE 

GO TU 110 
C 

c COMPUTE LINEAR SLOPES 

C 

90 00 100 nN=1 » NLS 

00 100 N3= 1 » 3 

COEP 1( NNt 3 f N3> =ALINE ( 2 * N31-AHNE 1 1»N3 ) 

COEP KNNt4,N3) = ALINE INN .N3) 

100 CONTINUE 
GO TO 130 
110 CONTINUE 
C 

C COMPUTE CUBIC SPLINE SLOPES 

C 

CALL SPFIT l MAXN *NLS » AL I NE »ELEN» COEP 1 » Kl» 0 » K2» EP »CP * K 3? K9) 

DO 120 N 3= 1 » 3 

COEP 1 1 NLS » 3 1 N3) *3. *COEP 1 1 NLS— 1,1 »N3 1 + 2 . ♦COEP 1 1 NLS- 1 , 2 , N3 1 +C0 EF II NL 
IS— it 3»N3 ) 

COEP 1 (NLS»4 ,N3) = C0EF l INLS-1 > 1 » N31 +C0EF1 (NLS- It 2 1 N3 1 +COEF l (NLS- it 3t 
l N31+C0EP 1 ( NLS— 1 1 4t N3 ) 

120 CONTINUE 
130 CONTINUE 

OJ 290 N=2 t NLD 

READ (101 KALlNE(NN*N3)tNN=ltNLSl tN3=lt3): 

IP INLS.LT. 3) GO TO 150 
NL=NLS— 1 
00 140 NN- 1 1 NL 

Ei=ABS(ALINE(NN, 1)— ALINE ( NN+lt 1 ) ) 

E2=ABS( ALINE (NN» 21 -AL I N£ ( NN+ l t 2 1 1 
E3=ABS(ALINE(NN, 31— ALINE(NN+1»3) ) 

IF (EUE2+E3.LE.EPS1 GO TO 150 
140 CONTINUE 
GO TO 170 

L50 DO 160 NN= 1 1 NLS 
DO 160 N3= 1 1 3 

C0EP2(NN,3iN31 = ALlNE(2t N3 l-AL 1NE ( 1 1 N31 
C0EP2(NNt4,N3J=ALlNE(NN,N3) 
loO CONTINUE 
GO TO 190 
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170 


CALC SPFIT IMAXN,NLS, ALINE, ELEN, C0EF2 ,K 1, 0 ,K2, EP , CP , K3, K4) 

00 180 N3=l,3 

CQEF21NLS,3,N3i=J.*C0£F2!NLS-l, 1 ,N3)*2.*COEF2I NLS- l, 2 , N3I+C0EF2 l NL 
1S-1, 3,N3J 

C0EF2! NLS,4,N3J=CUEF2INLS-l, 1 ,N3 ) *C0EF2 (NLS- 1, 2, N3i +C0EF2 t NL S~l , 3, 
IN3J+COEF2! NLS-1, 4,N3) 

180 CONTINUE 
C 

C STOKE PATCHES 

C 

190 DO 270 L=2,NLS 
DO 210 N.3= 1,3 
DU 200 M=l,2 
MM=MQO (M»2 ) 

LL=L-MM 

P ATCH I M , 1 , N31 =C0£F1!LL»4»N3) 

PATCHiM,2,N3)=C0EF21LL,4,N3) 

PATCHIM»3»N3]=SLUPE1 N-l ,LL,N3> 

PATCHIM,4,N3)=SLGPE! N,LL,N3) 

PATCH1M+2*1*N3J=C0EP11LL»3»N31 

PATCH!M+2,2,N3)=C0EF2!LL,3,N3) 

PATCH! M+2 » 3 » N3 ) = 0. 

PATCH! M+2»4,N3j=Q. 

200 CONTINUE 
210 CONTINUE 
C 

C COMPUTE PATCH IN FCKM CF S=MBM! T RANSPOSE I AND WRITE ON TAPE 

C 

DO 260 N 3= l * 3 
DO 230 14=1,4 
DO 230 J4=l , 4 
SUM=0. 

DO 220 K4=i,4 

220 SUM= SUM+XM AT i 14, -K 4 )*PATCHtK4»J4»N3) 

230 PAT! 14,J4)=SUM 
C 

DO 250 14=1,4 
DO 250 J 4= I , 4 
SUM= 0. 

DO 240 K4= 1,4 

240 SUM= SUM+PAT 1 14, K4) *XMAT I J4,K4i 
250 PATCH! 14, J4,N3)=SUM 
260 CONTINUE 

CALL RECUUT ! 7 ,2 , 0 ,P A TCH , 1 , 48 , 1 ) 

270 CONTINUE 
C 

C MUVE COEFFICIENTS 

C 

DO 280 N 3= 1 , 3 
DO 280 N4= 1,4 
DO 28Q NN= L , NLS 

COEF 1 1 NN , N4 , N3 1 = COEF 2 i NN , N4 , N3) 

280 CONTINUE 
290 CONTINUE 
RETURN 
END 
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Subroutine SPFIT 


Subroutine SPFIT uses a parametric cubic spline curve fit technique with optional 
enrichment of the given input curve. The method is explained in appendix A. The 
description, flow chart, and the FORTRAN statements for this subroutine are as follows: 

Language : FORTRAN 

Purpose : SPFIT is a parametric cubic spline curve fit subroutine. Parametric coeffi- 
cients are computed to approximate a cubic spline curve through a three-dimensional 

set of input points describing a curve, and, optionally an enriched curve is computed. 

Use : CALL SPFIT (MAXN, N, PNT, ELEN, COEF, NFIT, MAXSP, H, EPS, CPT, Kl, K2) 

MAXN The maximum number of input points allowed as stated in the 

dimension statement of the calling program. 

N The number of input points; 4 S N ^ MAXN. 

PNT A two-dimensional array of the consecutive points describing the 

three-dimensional (X,Y,Z) input curve. 

ELEN A one-dimensional array used by the subroutine for the chord lengths 

between each consecutive pair of input points. 

COEF A three-dimensional array used by the subroutine for the parametric 

cubic spline coefficients. 

NFIT A number of interpolated points to be computed between each pair of 

given points as specified by the user. 

MAXSP The maximum number of points allowed in the enriched curve as 

stated in the dimension statement of the calling program. If 
MAXSP is 0, only the cubic spline coefficients are computed and 
the calculation of the enriched curve is omitted. 

II The total number of points in the enriched curve calculated by the 

subroutine. 

EPS A small number supplied by the user which is used to check the sec- 

ond derivative at each point of the faired curve. The point will be 
omitted if the absolute value of the second derivative is less than 
EPS. An EPS of 0.0 will cause all the interpolated points to be 
retained. 

CPT A two-dimensional array used by the program for storage of the 

enriched curve. 
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K1 An integer supplied by the user. If Kl=l, retain all the input points. 

If Kl=2, include input points in second derivative test. 

K2 An integer supplied by SPFIT as an error code. If K2=l, normal 

return. If K2=2, error return when the number of interpolated 
points exceeds the allowable storage (MAXSP). 

Restrictions: SPFIT has been written with a variable dimension statement, and the 
following must be dimensioned in the calling program: PNT(MAXN,3), ELEN(MAXN), 
COEF(MAXN,4,3), CPT(MAXSP,3). If the coefficient-only option is used (MAXSP=0), 
dummy entries for NFIT, n, EPS, CPT, Kl, and K2 must be included in the calling 
sequence. The input curve must not include any consecutive duplicate points. 
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n r> o 


SUBROUT 1 Nt SPFIT ( MAXN, N , PisT , ELEN , COEF, NF I T , MAXS P , I i , EPS , CP T , Kl , K2 
1 ) 


C 

c 

c 

c 

c 

c 

c 

u 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


10 

c 

c 

c 

c 


20 


COMPUTES PARAMETRIC CUBIC SPLINE COEFFICIENTS TO 
APPROXIMATE A SMOOTH CURVE THROUGH A 30 SET OF INPUT 
POINTS ANO OPTIONALLY COMPUTES AN ENRl.CHEO CURVE 
MAXN IS THE MAXIMUM NUMBER OF INPUT POINTS ALLOWED 
N IS THE ACTUAL NO. CF INPUT POINTS 

NFIT IS THE NUMBER OF UESIRcO SPLINED POINTS BETWEEN GIVEN 
POINTS 

MAXSP IS THE MAXIMUM NUMBER OF SPLINED POINTS ALLOWED, 
MAXSP=(MAXN-l)*(MAX,NFIT + im FOR EPS OF 0. 

MAXS P= 0 OMITS COMPUTATION OF ENRICHED CURVE 
II IS THE NO. OF POINTS IN THE ENRICHED CURVE 
Kl IS AN INTEGER SUPPLIED BY THE USER 
Kl = l, RETAIN ALL INPUT POINTS 

Kl=2, INCLUDE INPUT POINTS IN SECOND DERIVATIVE TEST 
K2 IS AN INTEGER SUPPLIED BY SPFIT AS AN ERROR CODE 
K2=l, NORMAL RETURN 

K2=2, INCOMPLETE FAIRED CURVE WHEN MAXSP IS EXCeEOEO 

PROGRAMER - CHARLOTTE CRAIDON 2-1-71 

DIMENSION PNT l MAXN, 3 1 , ELEN I MAXN) , COEF ( MAXN ,4, 3) , CP U MAXSP , 3 > 

D1ST (X1,Y1,21,X2,Y2,22)=SQRT( ( X2-X1 > * *2+ 1 Y.2- Y1 )**2M Z2-Z 1 ) **2 ) 

N 1=N-1 

COMPUTE CHORD LENGTHS 
DO 10 NN=2 , N 

ELEN ( NN— l J =0 1ST! PNT ( NN— 1 , 1 1 , PNT I NN-i , 2» , PNT ( NN- 1 , 31 , PNTI NN , 1 1 ,PNT< 
INN, 2) , PNT INN ,3) ) 

CGNT I NUE 

SETUP COEFFICIENT MATRIX WITH UNCLAMPBD ENO POINTS 
(2ND DER=0. AT PI AND PNI 

COEFI 1,1,11=0. 

COEF I 1,1,21 = 2. 

COEF (1,1,3)=1. 

COEF ( N , 1 , 1 ) = 1 . 

COEF IN, 1,2 )=2. 

COEFIN, 1,31=0. 

DO 20 NN=2 , N 1 

COEF 1 NN , 1 , 1 ) =ELEN( NN> 

COEFINN, l,2)=2.*(ELEN(NN-li+ELEN(NN)) 

COEF (NN, 1, 3)=£LENINN-1) 

CONTINUE 


SOLVE FUR SLOPES 


DO 60 1=1,3 « 

COEFI 1,4, 11= I 3. /ELEN! U I *( PNT 1 2 , 1 J-PNT I 1 , 1 ) ) 

COEF ( N , 4 , 1 J = I 3 . / ELEN I N- 11 ) * ( P NT I N , I ) -PNT l N-l , 1 .1) 

DO 30 NN=2,N1 

COEF INN, 4, 1)=I3./(ELEN(NN-1)*ELENINN) ) > * l E LENI NN-1 ) **2 * { PNT ( NN+1 , I 
1 1-PNT INN , 1 1 1 +ELENINN 1**2* 1 PNT INN, I J-PNT 1NN-1-, 1 1 1 ) 

30 CONTINUE 
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c 

C SOLVE TRIOI AGONAL MATRIX 

C 

COEF ( 1,2 1 1 i = CUEE ( 1 , 1 1 3 J / CGEl- l it 1 »2I 
C0EF(lt3t li-CUEF ( 1, 4,1) /CLEF ( 1,1,21 
DO 40 K=2,N 
KM 1= K— 1 

TEMP=COEF(K, 1,21— COEFIK, It 1 i*COEF (KM 1 , 2 1 II. 

COEF IK 1 2 1 1 J-COEF (Kti»3)/IEMP 

CUEF (Kt3t 1 I— (COEFCKt 4t 1 1-CCEF(K, 1 , 1 J *COEFl KM 1 , 3, 1 1) /TEMP 
40 CONTINUE 

DO 50 K-= L t N 1 

KK=N-K 

COEF(KK,3t I J=COEF(KK,3, 1 1.-COEF ( KK ,2 , 1 1 *COEF l KK + i , 3 , 1 J 
50 CONTINUE 

ou CONTINUE 

C 

C COMPUTE CUBIC COEFFICIENTS FOR EACH SEGMENT 

C 

00 70 NN= 1 1 N l 
EL=1./ELENINN1 
EL2= EL*EL 
EL3=EL*EL2 
DO 70 1 = 1,3 

C0EF(NNt4t U =PNT INNt I ) 

E=PNT4NN*i,I I-PNTINN, I) 

C0EF(NN,2,l)=£*EL2*3.-EL*(2.*CGEF(NNt 3, 1 1 +.COEF (NN+l , 3 , 1 1 I 
COEF ( NN, 1 , 1 1 =— £*EL3*2. + EL2*( COEF (NN,3 1 1 l+COEF ( NN«-1 1 3 , 1 ) I 
C 

C REFERENCE LENGTH TO 1. 

C 

COEF INN t It I i*COEFI NN, It 1 1 /EL 3 
CUEFINN,2,I1=C0EF(NN,2,1I/EL2 
COEF ( NNt 3 , 1 i =COEF( NN ,3 , U /EL 
70 CONTINUE 

IF (MAXSP.EU.01 RETURN 
IF (Kl.EQ.01 Kl=l 
11=0 
C 

c compute enriched points 

c 

IF1T=NFIT+1 

XFIT=IFIT 

UELT=i./XFlT 

00 110 NN=1 , Ni 

00 100 NF = 1, IF I T 

E=NF-l 

T=0£LT*E 

IF (NN.EQ.l. ANO.NF.EW.il GO TO 80 
IF INF. EW.I. AND. Kl.EQ.il GO TO 80 
T6=6.*T 

EX=A85( T o*CUEF (NN, It 11*2. *C0EF(NN,2, 1 ) I 
EY=A8S(Tt>*C0EFlNN, 1 1 2 1 +2 . *C0EF ( NN 1 2 1 2 1 1 
EZ=A8Sl I&*CG£F(NN, 1, 3)«-2.*C0EF(NN,2,3) I 
EE-=(EX*-EY«-EZ)/lcLEN(NNl*ELEN(NN) I 
IF (EE.LT.EPSI GO TO 100 
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80 


90 

100 

110 

120 

130 


11=11*1 

IF 1 Il.GT.MAXSP) GO TO 130 
T2=T*T 
T 3=T *T 2 
00 90 1=1 t 3 

CPT1 II 1 1 )= T3*C0EF(NN ,1*1 1 *T2*C0EF CNN *2,11 *• T*COEF INN , 3* I ) *COEF (NN , A 

1 , 1 ) 

CONT INUE 
CONT 1N0E 
11 = 11*1 

IF C ll.GT.MAXSPJ GO TO 130 

00 120 1 = 1,3 

CPTC II »1 )=PNTIN» L) 

K2=l 
RETURN 
K2=2 
RE TORN 
END 
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Program ORTCON 

Program ORTCON (overlay (3,0)) is the control routine for the orthographic projec- 
tions of the input body. This program reads the plot information card and prints it, com- 
putes scale factors, computes vertical offsets for three-view plots, and notates on the 
plot. The flow chart and the FORTRAN statements for this overlay are as follows: 
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overlay (cbc, 3 , 0 ) 

PROGRAM UK i con 

c 

c CONTROL ROUTINE FOR ORTHOGRAPHIC PLOT'S 

C OF A SURFACE UR OF A COLLECTION OF SURFACES 

C 

COMMON /TnREEO/ ABODE IB) , hORZ.VERT .TEST l, PM I, THETA, PS I, 

1PLUTSZ.T YPE.NOU.NOrt, ISIuE ,KGOE 
C 

DIMENSION ORG ( 3 ) , ABC ( B ) 

DATA TYPED/ 3 HURT / , T YPEV/3FVU3/ 

C 

C 

C READ PLOT CARD 

C 

WRITE (6,10) 

10 FORMAT ( !H126X,27rtTriREE DIMENSIONAL PLOT DATA//) 

20 CONTINUE 

READ (5,30) ABODE 
30 FORMAT (BA10) 

IF (ENDFILc 5) 35, AO 

35 CALL NFRAME $ CALL CALPLT ( 0. ,0. ,999) $ STOP 

AO WRITE (6,50) ABODE 

50 FORMAT (1X.BA10/) 

DECODE l 72, 6 0, ABODE) HORZ , VERT ,T EST 1 , PH I , T HE T A , PS I , PLOTS/ , T Y PE , NOU 
1 , Nun , IS IDE , RUDE 

60 FORMAT ( 2 A2 , A3 , 3 F5 . 0 , 25 X , F5 . 0 , A3 , 3 I 3 , 7X , 1 1 ) 

IF (ISIOE.EQ.O) IS IDE=1 
C 

C READ PATCH TAPE 

c 

REWIND 7 

CALL RECIN (7,2,10, ABC, 1,8,1) 

IF (ENDFILE 7) 70,90 
70 WRITE (6,80) 

80 FORMAT (1H1/3BH ENO OF FILE ENCOUNTERED ON PATCH TAPE) 

STOP 

90 CONTINUE 

CALL RECIN ( 7 , 1 , IC ,XM IN , XMAX , YM1 N, YMAX, ZM1 N, ZMAX , NOB J ) 

IF (IS IDE* EU.2) YM IN=-Y MAX 
C 

C FIND SCALE FACTOR FROM MAXIMUM DIMENSION 

C 

XOI S=XMAX-XMIN 
YDIS=YMAX— YMIN 
ZD I S= ZMAX- ZM IN 
DMAX = AMAX1 (XD1 S.YDIS, ZD IS) 

SCALE=DMAX/PLGTSZ 
IF ( TYPE.NE. TYPE V) GO TO 1A0 
C 

C 3VU WHERE VIEWS ARE STACKED VERTICALLY 

C 

ORG( 1) =PHI 
ORG ( 2 ) =THET A 
ORG( 3)-=PSl 
PHI= THET A=PS 1=0. 
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YB IG=ORG ( 1 ) 

YORG-=FLOAT ( IFIX( YMAX/SCALE) ) *ORG( 1) 

IF ( YBIG.GT. GRG( 2) ) GO TO 100 
YB IG-ORG (2 ) 

YORG=FLOAT ( IF1X( ZMAX/SC ALE ) ) +0RGI2I 
100 IF (YBIG.GT. 0RGI3)) GO TO 110 
YBIG=0RG(3 ) 

YORG = FLO AT ( 1 FI X ( ZM AX /SCALE) ) +ORG (31 
110 CALL CALPLT (0.,YURG,-3I 
C 

C NOTATE ON 3VIEW PLOTS 

C 

NCHAR- IFIX(6.*PL0TSZ) 

IF (NCHAR. GT .BOJ GO TO 120 

X-0. 

GO TO 130 
120 CONTINUE 

NOI F= ( NCHAR— BO ) / 2 
X=FLOAT(NDIF)/o. 

NCHAR=BO 

130 CALL NOTATE ( X, 0 . , . 2 , ABC , C . , NCHAR ) 

YSAV=YMIN 

XMIN=YMIN=ZMIN=0. 

HORZ=lHX$VtRT=lHY 
YORG=URG( 1 )-YUKG-l 
CALL CALPLT (Q.,YURG,-3J 

CALL OTHPLT (XM1N, XMAX , YM IN , YMAX ,ZMI N ,Z MAX , NOB J , XM I 0 , YMID , ZM ID, SCA 
1LE) 

VERT = 1HZ 

YORG-ORG ( 2 )— ORG( 1) 

CALL CALPL1 (0. » YORG »-3 ) 

CALL JTHPL T (XMI N, XMAX , YM IN , YMAX , ZM IN , ZMAX , NOB J , XMIl) , YM ID , ZM ID , SCA 
1LE) 

HORZ= 1HY 

Y0RG=0RG(3)-0RG( 2) 

YM1i\I = (FlOAT( IFIX(YSAV/SCALE)-1) )*SCALE 
CALL CALPLT ( 0 . , YORG » — 3 ) 

CALL OTHPLT (XMI N, XMAX, YM IN , YMAX , ZM I N , Z MAX, NUB J , XM I 0 , YMI 0 , ZM ID , SCA 
lLt ) 

X=FL OAT (IF IX ( PLOTS Z+6.) ) 

Y=l. -0RG(3 ) 

GO T 0 160 
140 CUNT INUE 
C 

C CENTER PLOT 

C 

XMIU=. 5*(XMAX+XM IN) 

YM 10-= • 5*( YMAX* YM IN) 

ZM I 0= . 3* ( ZMAX+ZM IN ) 

XF I X = . 5* ( UMAX— XL) IS ) 

XM IN=XMIN-XF IX 
XMAX=XMAX+XF IX 
YF IX=.5*( DMA X- YD IS) 

YMIN-= YMIN-YF IX 
YMAX=YMAX+YF IX 
ZF1X=.3*(0MAX-Z01S) 

ZM1N=.ZMIN-ZF IX 
ZMAX = ZMAX+ZF I X 
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c 

c 

c 


c 

c 

c 


150 

C 

C 

c 


c 

c 

c 

160 


c 

c 

c 


ADJUST M1NIMUMS FOR GRID LINES 

XMIN=FLOAT(IFIX(XMIN/SCALE) )*SCALE 
Y Ml N= FLOAT (IFIXIYM IN/ SCALE! 1* SCALE 
ZMI W-=FLOAT ( I FI X( ZMIN/ SCALE ) ) * SC ALE 

NOTATE ID ON PLOT 

X=0. 

NCHAR=IFIX(li.*PLOTSZ)+3 
If (NCHAK.LE.80i GO TO 150 
N0IF= (NCHAR— 80) /2 
X=FLOAT ( NDIF ) / 1 1 • 

NCHAR=80 

CALL NOTATE ( X » • 8 t . 1 » A8C « C. » NCHAR ) 

CALL NOTATE (X» . 5 . . 1 » A8C0E ,0 . » NCHAR } 

ORTHOGRAPHIC 

CALL OTHPLT (XMlNt XMAX» YM IN» YMAX » ZMIN * ZMAX »NOB J »XM I D » YMI D» ZM ID* SCA 
1LE ) 

X= FLOAT ( 1 FIX ( PLOTSZ+2 . ) ) 

Y=0. 

END OF COMPLETE PLOT 
CONTINUE 

CALL CALPLT (X»Y»— 3J $ C A L L NFRAME 

IF (KODE.Ew.O) GO TO 20 

RETURN 

END OF GRTCUN 


END 
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Subroutine OTHPLT 


Subroutine OTHPLT determines the specified axis system and paper plane, sets up 
the rotation matrix, and establishes the necessary offsets for proper plot placement. The 
flow chart and the FORTRAN statements for this subroutine are as follows: 



SUBROUTINE OTHPLT IXMIN , XMAX , YM I N, YMAX , ZMI N, ZMAX , NOB J , XM I D , YM ID , ZM 
I I0,SCALtl 

ORTHOGRAPHIC PROJECTIONS 


COMMON/THREEO/ ABCOEI 8 ) » HORZ , VERT , TEST 1 , PHI , T HETA, P SI , 
l PLOT SZ t TYP E» NOU » NOW» I S I UE * KUOE 

OIMENS ION A12.3J ,NAME12J , ABCI8) 

DATA XSEE/2HX / t»YSEE/2HY /,ZSEE/2HZ /. 

1 XINT ST./3H0UT / » CONV/. 0 17 45 3293/ » NUN2/2 / » NAN2/ 24/ 

INITIALIZE 

ITEST1-1 

ITEST2=1 

If (XINTST.NE.TESTli 1TEST1=0 

IF ( PS I .cU.O. . ANU. THET A . EQ .0 * • ANU. PHI . EQ» 0* ) LTEST2 ; =0 

PHI=CUNV*PHI 

THET A=CONV *THETA 

PSI=CUNV*PSI 

SETUP AXIS 
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c 


S1NPSI=SINIPSI) 

SINTHE=SIN(THETA) 

SiNPHI^SINIPHII 
CUSPS I=CQSIPS1 ) 

COST HE-=COS (THETA) 

COSPHI=COS (PHI ) 

IF ( XSEE.NE.HURZi GO TO 20 
C 

C USE X FOR HORIZONTAL VARIABLE 

C 

IF ( 1 TEST2.EQ.G) GU TO 10 
A(1,1I=C0STHE*C0SPS1 

A 1 1,2)=-SINPSI*CUSPHI+S I NT HE*CGSPS I*S INPHI 
A(lt3)=SlNPS i*S 1NPH1 +SINTHE*COSPSI*COSPHI 
10 HMIN=XMIN 
HMAX=XMAX 
HMIU— XM10 
IHORZ-1 
GO TO 60 

20 IF ( YSEE.NE. HUKZi GO TO 40 

C 

C USE Y FUR HORIZONTAL VARIABLE 

C 

IF ( ITEST2.EQ.0) GO TO 30 
A( 1, l)=COSTHE*SlAiPSI ! 

A(1,2)=C0SPSI#CJSPHI+SINTHE*SINPSI*SINPHI 
A( It 3)=-C0SPSl*SINPHH-S 1NTHE*SINPSI*C0SPHI 
JO HMIN=YMIN 

HrtAX=YMAX 
HM10=YMI0 

IH0RZ=2 .j. 

GO TO 60 

c ! 

C USE Z FOR HORIZONTAL VARIABLE 

C } j ' 

40 CGNT INUE I 

IF ( I TEST 2. EU. 0) GU TO 50 
At it 1J=— SINTHE 
A( 1,2J=C0STHE*SINPH1 
Allf 3I=C0STHE*C0SPHI 
50 H.MIN = ZMIN 

HMAX=ZMAX 
HH1U=ZMI0 
IH0RZ=3 

60 IF ( XSEE.NE. VERT I GO TO 80 

C 

C USE X FOR VERTICAL VARIABLE 

C 

IF ( ITEST2.EQ.0) GO TO 70 
AI2. 1)=C0STHE*C0SPSI 
A(2t2i=-SINPSI*CQSPHI*-SINTHE*CGSPSI*S INPHI 
AI2t3i = SINPSI*SINPHI<-SI NTHE*COSPS l*COSPHl 
70 VMIN=XMIN 

VMAX=XMAX 
VMIU=XMIO 
I VERT= 1 
GO TO 120 

60 IF I YSEE.NE. VERT J GO TO 100 
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c 

c USE Y FOR VERTICAL VARIABLE 

C 

IF I I TEST2.E0.0J 00 TO 90 
At 2, i)=COSTHE*SINPSi 

At2,2J=CUSPSI*CUSPHH-SINThc*SINPSI*SlNPHl 
A( 2 1 3 ) -=— COSP SI*SINPHI+SINTHE*SINPSI *COSPHI. 

90 VM1N=YMIN 
VMAX=YMAX 
VM1D=YMI0 
1VEKT=2 

00 TO 120 
t 

C USE Z FOR VERTICAL VARIABLE 

C 

100 CONTINUE 

IF I ITEST2.Eu.0i GO TO 110 
At 2» II =-SI NT HE 
At 2, 2 ) =COS THE*S I NPHI 
At 2. 3i=C0STHE*C0SPHl 
HO VM IN— ZMI N 

VMAX=ZMAX 
VMIli=ZMID 

1 VERT- 3 

120 CUNT INUE 

C CENTER WITHIN PAGE SIZE IF SIZE GREATER THAN 28 INCHES 

C 

IF IPLUTSZ.0T.26..AND.TYPE.NE.3HVU3) VM IN=-1 3. *SCALE + FLO AT t I F I X t VM 
1 10/SCALE i ) *SCALE 
C 

C ROTATE MIDPOINT TO PLACE ROTATED VIEW CORRECTLY 

C 

IF t ITEST2.EU.0i 00 TO 130 

AM I D i= A ( It 1 J *XMI D+Al i,2)*YMI0+A(l,3i*ZMID 
AMID 2 = A t 2 » li *XMI D* At 2 » 2 ) *YMID*A(2»3 ) *ZMID 
HMIN-=HMIN— HMIO+AMl 01 
VMIN=VMIN-VM1D*AM102 
130 CONTINUE 
C 

c begin plotting 

c 

DO IbO I S I = I » ISIDE 
REWIND 7 

CALL RECIN l 7 t 2 * I C » ABC » 1 1 8 1 1 1 

CALL RECIN l 7 , I , IC t HI ,h 2 t H3 ,H4, H5, H 6 1 I 7 I 

DO 150 J = 1 » NOB J 

CALL RECIN 1 7 1 1 1 IC t N SURF » J3 » NAME 1 1) »NAME 1 2 ) » JA , J5 ) 

DO 140 N=i »NSURF 

CALL RECIN t 7 , 1 , IC »N01 , NSI , J3, J4, J5) 

CALL PLOT IT tNDl, NSI, I Si, I TE ST , 1 TESTl , IT EST2 , IHORZ.IVERT ,HMIN,VM1N 



I , SCALE , A J 

140 

CONTINUE 

150 

CONTINUE 

L60 

CONTINUE 


RETURN 


C 

C END OF OTHPLT 

C 

END 
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Subroutine PLOTIT 


Subroutine PLOTIT reads patch equations from tape and rotates them, computes 
enriched surfaces, and does a visibility test if desired. The flow chart and the FORTRAN 
statements for this subroutine are as follows: 
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SU OkuU T1 NE PLOTI T 1 NDl , NS1 » 1 S l » I TEST » I TEST 1 , I.TEST2 , 1 HORZ » I VERT « HMI 
IN, VMIN, SCALE, A) 

C 

C READS PATCHES FROM TAPE, 

C MANIPULATES IN SPECIFIED MANNER AND PLOTS 

C 

DIMtNS l UN PAT CHI 4, A, 31, PAT (4, 4, 21 ,A(2,3),PLPAT(4,21 , 

1PL1NE ( 34, 2i 

DIMENSION V£C(4»2J ,VPAT (41 

CUMMGN/THREEO/AbCUEI 81 , HORZ, VERT ,T EST 1 , PH I , T HET A , PS I , 

1PLOT SZ,TYPE,NOU,NUrf, ISIOE ,KUDE 
C 

NNU-NOU+2 

NNW=NUW+2 

FU-NUU+1 

FW=NOW+l 

DU= i ./FU 

DW=1 ./FW 

NPAT=NDi*NSi 

DU 230 N= 1 , NPAT 

CALC RECIN ( 7,2, 1C, PATCH, 1,48, 1) 

IF l 1 jI.E0.1J DU TU 20 
C 

C CHANCE Y SIGN 

C 

OU 10 14-1,4 
DO 10 J4= 1 » 4 

PATCH l I4,J4,2J=-PATCH( I 4, J4, 21 
10 CUNT INUE 

20 CONTINUE 

C 

C ROTATE PATCHES 

C 

IF l ITEST2.EQ.il GU TU 40 
UU 30 14=1,4 
DU 30 J4=l ,4 

PAT ( 14, J4, 1 1 =PAT CH 1 14, J4, IHuRZl 
PAT ( 14,J4,2)=PAT CH( 14,J4,1VERT) 

30 CONTINUE 

GU TO 80 
40 CONTINUE 

DO 70 14=1,4 
DO 70 J4=l ,4 
DO 60 K2=l , 2 
PAT l 14, J4,K2 1=0. 

DU 50 N3=l , 3 

PAT ( 14, J4,Kt J = PAT( 14 , J4 , K2 1 + A (K2 ,N3 I * PATCH 1 I 4, J4, N3J 
30 CONTINUE 

60 CONTINUE 

70 CONTINUE 

80 CONTINUE 

C 

C PLOT IN W DIRECTION 

C 

DO 150 NU= 1 , NNU 
EU=NU— 1 
U=EU*DU 
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DO 90 J4= 1,4 
00 90 K2= 1 » 2 

PLPAT ( J4»K.2) = ( (U*PAT 1 1 , J4,K2 ) +PAT42, J4,K.2fc)*U+PAT( 3, J4»K2 ) )*U+PAT( 
14, J4,K2> 

VEC4 J4,K2) = 4 3.*U*PAT l 1 , J4,K2 ) +2 . *PAT ( 2 , 04, K2 ) )*U+PAT (3* J4,N2 ) 

90 CONTINUE 
NIT=0 

00 140 Nw = 1 » NNW 
E W=NW— 1 
lni=Ert*0W 

IF ( ITESTl.EU.O) GO TO 12C 
C 

C COMPUTE DV/OU AND DV/Dw 


C 

DO 100 0=1,2 

VP AT (0) = U W*V EC ( 1, 0) +VECI 2,0 ) ) *VEC ( 3 , 0 ) i *irf*VEC ( 4, 0 ) 

VPAU J*2J = 4 3.*PLPAT( 1 , J )*N+2.*PLPAT (2,01) ♦•ti+PLP AT 13,0) 

100 CONTINUE 

VNuRM= VP AT (i)*VPAT44)-VPAT(2)*VPAT(3) 
v . "~lf H-Sl.Ea.l) VNORM=-VNQRM 
IF 4VNURM.GE.0.) GO TO 120 
IF (NIT.GT.l) GU TO 110 
N l T = 0 
GO TO 140 

110 PLINE4N1T«-1,1)=HMIN$PLINE4N1T*1,2)=VMIN 
PLINE(NIT+2,1)=PLINE(N1T*2,2)=SCALE 
CALL LINE 4PL1NE41,1) , P L I N E l 1,2) , NIT, 1,0, 0,0) 

Nl T=0 
GO TO 140 
120 NiT=NLT+l 

DO 130 K.2=i,2 

PL1NEINIT,K2) = ( I W*PLPAT ( i,K2)+PLPAT 12 , K2 ) J *W +PLPAT I 3 ,K2 ) )*W+PLPAT( 
14.K.2) 

130 CONTINUE 
140 CONTINUE 

IF (NIT.LE.l) GU TO 150 

PLINEI NlT+1, 1)=HMIN*PLINE(NIT+1,2)=VMIN 

PLlNE (NIT +2 , 1 )=PL1NE( NI T +2,2)=SCALE 

CALL LINE ( PLINE ( 1, 1 ) ,PL INE( 1,2) ,NIT , 1 ,0,0,0) 

150 CONTINUE 
C 

C PLOT IN U DIRECTION 

C 

DO 220 NW= 1 , NNW 
Ew=N*-l 
is = EW*DW 
DO 160 J4= 1 , 4 
DO 160 K2= 1 , 2 

PLPAT (J4,K2)=( 4**PAT4J4, 1,K2)+PAT( J4,2,K2) )*W+PAT( J4,3,K2) )*W+PAT( 
1 J4,4,K2) 

VEC4 J4,K2)=(3.*w*PAT4 J4, t,K2)*2.*PAT( J4,2-,X2))*W*PAT4 J4,3,K2) 

160 CONTINUE 
NI T=0 

DO 210 NU= 1 , NNU 
EU=NU— 1 
U=£0*DU 

IF 4 ITESTi.EQ.O) GO TO 190 
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c 

C COMPUTE OV/DU AND DV/DW 

C 

DO 170 J=i,2 

VP AT 4 J) = < 3.*PLPAT( 1, J ) * U+2.*PLPAT (2 , J) 1*U*PLPAT(3, J) 

VP AT ( J +2 J - ( ( U*V EC ( 1,J)+VEC(2,J) ) *U+VEC ( 3, J i ) *U«* VEC C 4 * J ) 

170 CONTINUE 

VNORM=VPAT m*VPAT(4)-VPAT(2)*VPATC 3) 

if (ISI.EO.il VNCRM=— VNCRM 

If l VNORM.GE.O.J GO TO 190 

If tNIT.GT.il GO TO 180 

N1 T= 0 

GO TU 210 

180 PLlNEt NiT+1 , i)-HMIN$PLINE(NIT+l,21=VMlN 

PLINEtNI T + 2, il=PUNElNi T*2,21 = SCAlE 
CALL LINE ( PLI NE ( 1 » 1 ) t PL INE ( L#2i »N IT 1 1 f Ot UifO ) 

N l T= 0 
GO TO 210 
190 " N I T-N IT + i 

DO 200 K2= 1 f 2 

PLINEtNIT ,K2)=((U*PLPAT ( 1,K2J+PLPAT12»K2)1*U+PLPAT( 3 , K2 ) 1 *U+PLP AT ( 
14, K2) 

200 CONTINUE 
210 CONTINUE 

If CNIT.LE.il GO TO 220 

PLINECNIT+l, 1 ) =HM N* PL I NE(NIT*1, 2)^VMIN • 

PL INE iNl T + 2 V U-PLINEC NIT + 2» 21-SCALE 
CALL LINE CPLlNECLvllvPLINEC 1,2) ,N1T , 1 , 0, 0, 0 1 
220 CONTINUE 
230 CONTINUE 
RETURN 
END 
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Program XPLT 


Program XPLT (overlay (4,0)) is the control routine for cross-section plots through 
the input body. The program reads the plot information card and prints it and notates on 
the plot. The flow chart and the FORTRAN statements for this overlay are as follows: 
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OVERLAY (CBC, 4,0) 

PROGRAM XPL T 
C 

C CONTROL PROGRAM FOR CROSS 

C SECTIONAL OR CONTOUR PLOTS 

C 

CUMMGN/XSECT/ABCUEI6) ,PPL 1(3),PPL2(3) ,PPL3(3) , 
IPLGTSZ,HPAGE,VPAGE, INP,NCU,NOW, 1 SIDE, I PR IN. 
2K.00E.XSTAT .THETR.XMACH 
DIMENSION ABC ( 8 J . ABCDI 8 I 
C 


c 

c 

read plot cardis) and print 

WRITE (6,10) 


10 

FORMAT (1H127X, 25HCR0SS SECTIONAL PLOT 

DATA//) 

zo 

CONTINUE 

READ (3,30) ABCUE 


30 

FORMAT (8A10) 

IF (ENOFILE 5) 35,40 


35 

CALL NFRAME $ CALL C ALPL T ( 0 . ,0. , 999 ) $ 

STOP 

40 

CONTINUE 
WRITE (6,50) 


50 

FORMAT (/26X,30H ******* PLOT CAROIS) 

*«*****// ) 


WRITE (6,60) AbCDE 


60 

FORMAT ( IX 8 A 10 ) 



DECODE (80, 70, ABCDE) X1,Y1,Z1,X2,Y2,Z2, 
1 , 1 NP , NUU , NOW , I CU T , I S I DE , IPR 1 N , RODE 

X3 t Y3 tZ3# PLOTSZ tHPAGfc, VP AGE 

70 

FORMAT (iOF6.0,2F3.0,A3,2I3,I2,3Il) 
IF ( ICUT . EG. 0) GO TO 90 



READ (5,30) A8CD 
WRITE (6,6 0) At) CD 
DECODE (21 ,60, ABCD) DX,DY,DZ,IH 
80 FORMAT (3F6.0.I3) 

90 CONTINUE 

IF ( INP.EG.3HANG) GO TO 110 

WRITE (6,100) X1»Y1, Z1,X2 »Y 2, Z2 » X3 »Y 3 »Z 3 .PLUTSZ, HP AG £»VP AGE 
100 FORMAT (//13X.1JHCUTTING PLANE/ 6X , 1HX l 1 X, 1HY 11 X, 1HZ/3F 12. 5/3F12. 5/ 
13F12.5//6X,5HSCALE7X,5HHPAGE7X,5HVPAGE/F12.5,2F12.2) 

GO TO 130 

110 WRITE (6,120) XI , Y 1 , Z1 , X2 , Y2 , Z2 , PLOTSZ, HPAGE , VP AGE 
120 FORMAT (// 13X , 13HCUT TING PLANE/6X,2HX010X, 2HY0 10X, 2HZ0/3F 12. 5/6X.1 
IHX11X,3HTHETA7X,4HMACH/3F12.5//6X,5HSCALE7X,5HHPAGE7X,5HVPAGE/F12. 
25.2F12.2) 

130 CONTINUE 

WRITE 16,140) NUU, NOW, ICUT , IS IDE , IPRI N 
140 FORMAT (6X,28HN0U NOW 1C0T LSIOE I PR IN/1.9 , I 5 , 1 6, 2 1 7) 

IF ( IS IDE. EQ • 0 ) ISI0E=1 
IF ( IPRI N. EQ . 0 ) IPRIN=1 
HSA V=HPAGE $ I F ( IH .EQ. 0 ) HP AGE=0. 

NCUT=ICUT+1 

IF ( ICUT.NE.O) GO TO 130 
DX=0.$DY=0.$DZ=0. 

GO TO 200 

150 *RITE (6,160) DX,DY,DZ 

160 FORMAT l 6X , 2 HDX1 OX ,2HDY 1 OX , 2HDZ/ 3F12. 5) 

IF UH.NE.O) GO TO 180 
WRITE (6,170) 

170 FORMAT ( 27H OVERLAID PLOTS WITH IH = 0) 

GO TO 200 
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o o o 


180 

190 


200 


210 

220 

230 

290 

230 

2t>0 

270 

280 

290 


C 

C 

c 


300 

310 

320 

C 

C 

C 


330 


WRITE (o,l90) 1H 

FORMAT (29H SPACED PLOTS WITH IH * ,13) 
LOUP FOR INCREMENTED CUTS 


00 350 N=1»NCUT 
IF IN.EU.l) GO TO 290 
IF ( INP.EU.3HPNT ) GO TO 250 
a2=X2+OX$Y2=Y2+OY$22=22+OZ 
DECODE ( 10,210, A8C0E (2) ITEM PI 
( A8 ) 

110, 220, ABODE (9) 

16XA9) 

(90»^30»A8CDE(2) 

( A6 »3Fo.2 » A9) 

(6,290) A2 , Y2» 22 


) TEMP2 

)TEMPi,X2,Y2,22,TEMP2 


CUTT ING PLANE/ 6X» 1HX11X, 5HTHETA7X* 4HMA 


FORMAT 
DECODE 
FORMAT 
ENCODE 
FORMAT 
WRITE 

FORMAT (// ZX»25H INCREMENTED 
1CH/3F12.3) 

GO TO 290 

Xi=XH-OXiYL=YU-DY*Zl=Zl + DZ 
X2=X2-kDXS.Y2=Y2+DY$Z2 = Z2+0Z 
X3=X3*UXFY3= Y3+U Y3 2 3=23+02 
DECODE ( 10 ,260, A8C0E16) JTEMP1 

ENCODE ( oO» 2 70 , AUCDE ( 1 ) )X 1 , Y l , Z 1 »X2 , Y2 , 22 , X3* Y3 , Z3, T tMP 1 
FORMAT ( 9F6. 2, A6) 

WRITE 16,280) XI , Y1 , Zi , X2 , Y2 , Z2 , X3 , Y3 , Z3 

FORMAT ( // 26H INCREMENTED CUTTING PL ANE /6X, 1HX 1 IX » IH Y 1 IX, 1HZ/ 3F l 2. 
15/3F12.5/3F12.3) 

CONTINUE 

PPLi(l) = Xl$PPLl(2)-YUPPLl(3)=Zi 
PPL2 ( 1)=X2*PPL2( 2)=Y2$PPL2(3)=Z2 
PPL3 ll)=X3$PPL3(2)=Y3J>PPL3(3)=Z3 
IF (N.EQ.NCUT) HPAGE=HS A V 


READ PATCH TAPE 


REWIND 7 

CALL RECIN ( /, 2 , I C , ABC , 1 , 8 , 1 ) 

IF (ENDFILE 7) 300,320 

FORMAT (IH1/38H END OF FILE ENCOUNTERED ON PATCH TAPE) 
STOP 

CONT INUE 

NOTATE 


X-0. 

IF ( HPAGt.Ew.O. ) GO TO 330 
NCHAR=IFIX( ll.*HPAG£H-2 
IF (NCHAR.GT.80) NCHAR=80 
CALL NOTATE ( X , . 8 , . I , ABC , 0 . , NCH AR ) 

IF UCUT.NE.O) CALL NOTATE (X,.6,.1,ABCD,0.,NCHAR) 
CALL NOTATE ( X , . 9, . 1 , ABCDE , 0. , NCHAR ) 

CONTINUE 

Y=FLOAT( IFIXl .5* VP AGE) )+l. 

CALL CALPLT (X,Y,-3) 

I F ( HPAGE.NE.O. ) CALL NOTATE (0. » 0. » 2. » 3* 0 • »—l ) 
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CALL xcur 

x=o. 

y--y 

CALL CALPLT (0.»0.»3) 

CALL CALPLT iX,Y,-3J 
if IHPAGE.EO.O.i GO TO 340 
X=HP AGE+2. 

Y=0. 

CALL CALPLT (X,Y t -3) 

CALL NFRAME 
340 CONTINUE 
350 CONTINUE 

IF iK.OOE.Ew.Oi GO TO 20 
RETURN 
C 

L END OF XPLT 

C 

ENO 
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Subroutine XCUT 


Subroutine XCUT sets up the transformation matrix and the origin of the plot 
coordinate system for the cross-section plots. The flow chart and the FORTRAN state 
ments for this subroutine are as follows: 



subroutine xcut 

CROSS SECTIONAL OR CONTOUR PLOTS 

COMMON/X SECT/ ABC OE (8J#PPL1(3I iPPL2(31fPPL3(3J» 

1 PLOTS/ »HPAGE» VPAGE» I NP » NCU t NOrt j I SIDfc, I PR IN * 
2KU0E.X3TAT .THETK.XMACH 
DIMENSION AI 3» 3) » NAMEt 2 ) t ABC ( 8) 

DIMENSION XPRM(J),PNTH31*PNT2(3),PNT3{3),PNT4(3J 
EQUIVALENCE IPNT i» PPL1 J , IPNT2,PPL2) , (PNT3tPPL3 J 

IN IT l AL 1 2c 

SCALE=1. /PLOTS/ 

Xd=PPLH L) »YB=PPL1C2)$/B=PPLH3) 

IF l INP. EQ.3HANGJ CO TO 20 
M=lPCOFl ACOEFt B»C»D,PPL 1 ,PPl2,PPL3) 

IF (M.EQ.i) GO TO 30 






n C' r o r r o o n 


u 


WRITE ( 6 i 1 Oi 

FORMAT (// 19H NO PLANE OESCR IBEO// ) 

RETURN 
20 XSTA T=PPL2 ( 1 J 

THET R=PPL2 (2)*. 01745329252 
XMACH-PPL2 ( 3 ) 

Y^-Y 

THE TR=PPL2<2i *.01745329252 
XMACH=PPL213) 

BETA^SURT ( XMACH**2-1 . ) 

ACOEF-i. 

B=-BETA*COS( THETR) 

C=-BETA*SINl THETR) 

D=XS TAT 
PNTi ( 1 )=XSTAT 
PNT 1 I 21 - 0 . 

PNTI (3 ) s 0. 

PNT2(2)=PNT3(3)=G. 

PNT2(3)=PNT3I2i=l. 

PNT2( L)=PNTH 1)-C 
PNT3 C l )-PNT 1 ( 1 ) — B 
30 CONTINUE 

PNT4I l.)=PNT3(l)4>PNT4(2)=PNT3(2)$PNT4(3)=PNT3(3) 

COMPUTE VECTORS 

T IX-PNT3I l )-PNT l ( l ) $ T1Y=PNT3 (2 1— PNTI (2 ) ST l.Z= PNT3 (3 )— PNT 1 1 3 ) 
T2X=PNT4(i)-PNT2(l)$T2Y=PNT4(2)-PNT2(2) ST2 2= PN T4 ( 3)-PNT2 ( 3) 
FNX=T2Y*T1Z-T1Y*T2Z 
FNY=IiX*T22-T2X*TLZ 
FNZ=T2X*T1Y-T1X*T2Y 
UN=SURTIFNX*FNX+FNY*FNY+FI\Z*FNZ) 

UNX=FNX/UN 
UNY=FNY/UN 
UNZ-FNZ/UN 

UT 1= S0R1 IT1X*TIX+T1Y*TIY+T1Z*I1Z) 

UT IX=T1X/UT i 
0TlY=TilY/UTI 
UTiZ^TIZ/UTl 
uT2X=UNY*UT1Z-UnZ*UTIY 
UT2Y=UNZ*UTIX-UNX*UT1Z 
U T2Z =UNX*U T 1 Y- UN Y*UTiX 

SETUP TRANSFORMATION MATRIX 

A( 1. 1)=UNX$AU,2)=UNY$AI l,3)=UNZ 
AI2tli-UTlX$A(2t2)=UTlY$A(2f3)=UTlZ 
A(3t i)=UT2X$A(3,2)=UT2Y$AI 3,3)=UT2Z 

SET ORIGIN OF NEW COORD. SYSTEM 

XPRMl 1 )=Al 1, 1)*XG+AI i» 2i *YB+AT1 »3)*ZB 
XPRH 12 ) = A( 2 » 1)*XB+A12»2)*YB+A(2»3)*ZB 
XPRM(3)=A(3»l)*XB*A{3»2)*Y3+A(3,3)*ZB 
HMIN=0. 

VMIN=0. 

WRITE (6,40) ACOEFfB.CfD 
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40 FORMAT (33H EQUATION OF THE PLANE AX+6Y+CZ=0/3H A=El 5. 8 ,4X ,2HB= » EL 
I5.a,4X,2HC=,E15.8,4X,2H0=»E13.8) 

C 

C bEGIN COMPUTING ANO PLOTTING LINES OF INTERSECTION 

C 

00 80 I SI- l » I S IOE 
REMIND 7 

CALL RECIN l 7, 2, 1C , ABC , l » 8, 1 J 

CALL RECIN ( 1 , 1 , IC,Hl,H2,F3,H4,H5,H6,N0BJI 

00 70 J = I » N03J 

CALL RECIN ( 7# It 1C »NSURF » J J » NAME ( 1 1 1 NAME! 2 ) » J J » J J ) 

WRITE (6,50) NAME 
30 FORMAT (IX2A10) 

00 60 N= UNSURF 

CALL RECIN ( 7, U IC,N01,NSUJ3, J4, J5) 

CALL PLANEX CND1 »NSI* IS 1 , ACOEF , b,C » D, HM IN, VMIN , A , SC ALE , XPRM, NAME ) 
60 CUNT INUc. 

70 CONTINUE 

80 CONTINUE 

RETURN 
C 

C END OF XCUT 

C 

END 
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Subroutine PLANEX 


Subroutine PLANEX reads patch equations from tape and finds the patch equation of 
the line of intersection of a plane and the patch. A line of intersection through the patch 
is computed if such a line exists. The flow chart and the FORTRAN statements for this 
subroutine are as follows: 



54 









n r n n o 


SUBROUTINE PLANEX ( ND 1 , NSl , l S l , ACOEF, B , C.O.HMI N, VMI N, A, SCALE , XPRM, 
1 NAME ) 


READS PATCHES FROM TAPE, 

SOLVES FOR INTERSECTIONS BY SPECIFIED PLANE 
AND PLOTS 


C 


c 


c 

c 

c 


10 

20 

c 

c 

c 


so 

c 

c 

c 


40 


so 


DIMENSION PATCH! 4, 4, 31 , A ( 3 , 3 1 , AL I NE ( 52 , 3) , 

1GMAT (4,4) » GV EC l 4 ) »VEC(4»3) 

DIMENSION XPRM(3) ,NAME1 2) 

COMMON/XSECT / ABODE ( 8 ) »PPL1(3)»PPL2L3) ,PPL3(3i, 

1 PLOT SZ» HP AGE , VP AGE, I NP ,NGU,NOW, I SLOE, IPRIN , 

2KOOE.XSTAT ,THETR,XMACH 

NNU=N0U+2 

NNW=N0w+2 

FU=N0U+1 

FW=N0W+1 

du=i./fu 

0W=1 . / FW 

NPAT=ND1*NS1 

DU 370 NNPAT =1 ,NPAT 

CALL RECIN ( 7,2, IC, PATCH, 1,48, II 

IF l ISI.EQ.l) GO TO 20 

CHANGE Y SIGN 

DU 10 14=1,4 
DO 10 J4= 1 ,4 

PATCH l 14, J4, 2)=-PATCH( 14, J4, 21 

CONTINUE 

CONTINUE 

SOLVE FOR G MATRIX 

DO 30 N4=l ,4 
DO 30 14=1,4 
GMAT (N4, 14J=0. 

GMAT ( N4 , 14)=ACGEF*PATCHIN4, 14,1 ) + 6*PATCH(Ni4, 14 ,2 ) +C*PATCH( N4 , 14,3) 
CUNT I NOE 

SOLVE FOR W WITH U=0. 

U=0. 

DO 40 1=1,4 
GVEC 1 1)=GMAT (4, I ) 

GVEC 14)=GVEC14)-D 
GUESS=0. 

IKUDE=KUdSOL l GUESS ,W ,GVEC) 

GO TO (SO, ISO, 150) , IKODE 
NPT= 1 

CALL VSOLV (X,Y, Z,U,W, PATCH) 

ALINE ( NPT , 1 ) = XS>AHNE (NPT ,2) = YiAL INEINPT ,3)=Z 
DO 140 N=2 , NNU 
E=N- 1 
U=E#0U 
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c 

C SOLVE FUR cubic IN « 

C 

DU 60 J4=l < 4 

GVEC ( J4 ) = l IU*GMAT( 1, J4) + GMAT (2»J4))*U«-GMAT(3»J4))*U+GMAT(4,J4) 
60 CONTINUE 

GVEC(4)=GVEC(4)-0 

GUESS=W 

IKODE=KUBSOL ( GUESS ,*, GVEC) 

GU TU (' 1 30 » 70 » 70) » IKUuE 
C 

C TRY U EUR *=i. ON OUTER 30UNDARY 

C 

70 w=i. 

DO 90 1=1,4 
GVEC ( 1 )=0. 

00 bu J= 1 , 4 

60 GVEC 1 1 ) =GV EC ( I ) + GM AT ( 1 , J ) 

90 CONTINUE 

GVECI4I=GVEC{4J-U 

GUESS=i. 

I KOUE=KUBSGL (GU ESS, U, GVEC) 

GO TU (120,100,100), IKOOE 
C 

C TRY U FOR W=0. ON OUTER BOUNDARY 

C 

100 W = 0. 

00 110 1=1,4 
110 GVEC ( I ) =GMAT (1,4) 

GVEC ( 4 ) =GV EC ( 4) -0 
GUtSS=i. 

IKOOE=KUBSOL (GUESS, U, GVEC) 

GU TO ( 120 , 3o0 , 360 ) , IKODE 
120 CALL VSOLV ( X , Y , Z, U, M , PATCH) 

NPT=NPT+1 

ALINE (NPT, 1)=X*ALINE(NPT,2)=Y$ALINE(NPT,3)=Z 
GO TO 3o0 

130 CALL. VSOLV IX»Y»Z»U»W, PATCH) 

NPT=NPT+i 

AL1 N t(NPT»l)=XSALINE (NPT, 2)=Y SALINE (NPT* 3)-=Z 
140 CONTINUE 
GU TO 360 
C 

C SOLVE FOR U WITH W=0. 

C 

150 W=0. 

DO 160 1=1,4 
160 GVEC(l)=GMAT(I,4) 

GVEC (4)=GVEC(4)-D 
GUESS=0. 

I KQ0E=KU3SUL ( GUESS ,U,GVEC) 

GO TO 1 170,250,250), IKODE 
170 NPT=1 

CALL VSOLV (X,Y»Z»U,W»PATCH) 

ALINE (NPT, 1)=XSALINE(NP T,2) = Y SALINE (NPT, 3):=Z 

DO 240 N=2 «NNW 

E=N-1 

W=t*DW 
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j o u ! -> o o 


SOLVE FOK CUBIC IN U 


DO 130 J4= 1 , 4 

GVEC I J4) = I I W*GMAT I J4 • 1 ) +GMAT ( J4 « 2 ) ) *W +GMM I J4 , 3) ) *W+GMAT I J4 » 4 ) 
CUNT INUE 

GVEC (4)=GVEC 14)-U 
GUE3S-U 

IKOOE = KUaSGL(GUESS»U »GV£C) 

GO TO (230,190,190), IKODE 

ray w for u=i. on outer bounoary 


190 0=1. 

UO 210 I = L , 4 

GVEC ( l )=0. 

DO 200 J = 1 ,4 

200 GVEC ( I )=oVEC 1 l) + GMAT I J, I » 

210 CONTINUE 

GVEC(4)=GVEC(4)-D 

GUESS=1. 

1KUDE=KUBSUL ( GUESS, WjGVECI 
GO TO (220,360,3601, IKODE 
220 CALL VSOLV ( X ,Y , Z,U, W,P ATCH) 

NPT=NPT+1 

AL INE l NPT* 1 ) = XSAL 1 NE I NP T , 2) = Y$AL INt ( NPT , 3 )-=2 
GO TO 360 

230 CALL VSOLV I X ,Y , 2,0, W, PATCH) 

NP T=NPT+ 1 

ALINE (NPT, 1)=X$ALINE(NPT ,2 l = V $AL INE ( NPT ,3 i-Z 
240 CONTINUE 
GO TO 360 
C 

C SOLVE FOR U WITH W= l . 

C 

250 W=l. 

DO 270 1=1,4 
GVEC ( 1 ) =0. 

DO 260 J=1 ,4 

260 GV£C( l)=GVEC(i)«-GMAT(I» Ji 
270 CONTINUE 

GVEC (4)=GVEC (4)-D 
GUESS=0. 

IKGDE=KUBSUL(&UESS»U*GVEC) 

GO TO (280,370,370), IKODE 
2B0 NPT=1 

CALL VSOLV (XfYfZfU, W, PATCH ) 

ALINE! NPT, 1) =X$ ALINE (NPT,2)=Y$AL INE (NPT ,3 ) -L 
00 350 N=2, NNW 
E = N— 1 
W= l.-E*Dw 

SOLVE FOR CUBIC IN U 


DO 290 J4=i » 4 

GVEC ( J4) = I ( W*GMA T ( J4 , 1 ) +GMAT I J4,2 


♦W + GMAU J4,3) )*W+GMAT( J4,4) 


C 

C 

c 



290 CONTINUE 

GVEC 14) =G V EC ( 4 J - 0 
GUtSS=U 

lKQDc=KUB30L( GUESS, U, GVEC) 

GU TO ( 340 ,300,300) , IKOOE 

C 

c TRY H FUR U-l« ON OUTER BOUNDARY 

C 

300 U- 1 • 

DU 320 1=1,4 
GVEC( U = 0. 

DO 3 10 J=l,4 

310 GVEC ( I )=GVEC ( I ) + GMAT i J , U 
320 CONTINUE 

GV£C(4)=GVEC(4)-U 

GUESS^l. 

IKGUE-KUBSGL (GuESS»W,GVEC) 

GU TO ( 330 y 360 , 3 60 J , IKGUE 
33U CALL V SOL V ( X , Y , Z , U , W , P A TCH) 

NPT=NPT+1 

ALINE (NPT, 1 ) =X$ALI NEC NP T, 2 )=Y SALINE (NPT ,3L=Z 
GO TO 360 

340 CALL V SOL V ( X , Y f Z t U , w , P A TCH ) 

NP T-= NP T + 1 

ALINE (NPT, li^X$ALINE(NP f , 2)=Y$AL INE ( NPT , 3) =Z 
330 CONTINUE 
C 

C SCALE, PRINT AND PLOT 

C 

3o0 CUNT I NUE 

IF (NPT.LE.I) GU TO 3/0 

CALL PLTROT (NPT , A , AL INE , HM IN, VMI N , SC ALE, ; IPR IN ,XPRM, NAME) 
370 CONTINUE 
RETURN 
C 

C END OF PLANtX 

C 

END 
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Subroutine VSOLV 


Subroutine VSOLV evaluates a patch equation for x-, y-, and z-coordinates with u 
and w given. The flow chart and the FORTRAN statements for this subroutine are as 
follows: 


( Subroutine VSOLV 


5 

7 

Solves 

|yj = V(u,w) for 
u and w given 


Z 


SUdRUUT I Nt VSOLV t X , Y , l , U , W , PATCH) 

C 

C SOLVES FOR VCU.W) FROM PATCH EQUATION 

C 

DIMENSION VEC(4,3J , PATCH 1 4 , 4 , 31 ,V(3J 
C 

DU 10 J 4 = 1 1 4 
DU 10 R3=l»3 

VECl J4,K3)=( (U*PATCH(1, J4 , K 3 J *PA TCH 1 2 , J4, K3J ) *U+PATC HI 3 , J4, K 3) ) *U«- 
1PATCHC4, J4,K3) 

10 CONTINUE 

DU 2 0 

VU3)=UW*VECllt K3J+VECt2tK3JJ*MtVEC( 3, K3j J*W*VEC ( 4, K3) 

20 CONTINUE 

X=Vt 1)$1T-VI2I*2=V(3) 

RETURN 

END 
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o o o 


Function KUBSOL 


Function KUBSOL selects the required real root from the roots of a cubic equation. 
The flow chart and the FORTRAN statements for this function are as follows: 



FUNCTION KUBSuLl IL,T .CUtFFSi 
C 

C FINOS THE ROOTS OF A CUBIC AND SELECT Si THE REQUIRED REAL 

C ROOT BETWEEN 0. AND 1. CLOSEST TO A GIVEN ESTIMATE 

C 

C THE ROUTINE SETS 

C KUB SGL= 1 SUCCESS 

C KU6S0L= 2 NO ACCEPTABLE ROOT 

0 KU6SUL= 3 ERROR 

C Tl IS ESTIMATE FOR SELECTION OF REQUIRED ROOT IF ALL ARE REAL 

C AND BETWEEN 0. AND i. 

C T IS REQUIRED ROOT WITH KUBS0L=1 

C CDEFFS ( 1 ) AkE THE COEFFICIENTS AT**3+BT**2+CT«-0=0. 

C 

DIMENSION CUEFF3U) 

COMPLEX R00TSI3) ,FTEM(8) 

DATA EPS/i.t-6/,EPi/l.E-5/ 

CHECK FOR GIVEN VARIABLE ON SEGMENT END POINTS 


T = 0. 

IF (ABSIC0EFFSIAJI.LE.EPS) GO TO 90 
T = l. 

CCC=C0EFFS(i) + CUEFFS(2)*CGEFFS 1 3) +CUEFF S (4-) 
IF IABSICCCJ.LE.EPS) GU TO 90 
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o n 


c 

c 

c 

10 

20 

30 

40 


SOLVE CuiJIC EQUATION FOR T 

IF (ABS(CGEFFS( 1)) .GT .£P 1 ) GCJ TU 40 
IF l ABS(CGEFF3(2 ) ) .GT.EP1 ) GO TO 30 
IF I ABSl CGEFFSl 3) I .GT .EP 1 ) GO TO 20 
KU6SUL=3 
RETURN 

T=-Cu£FF3 ( 4) /CGEFFSl 3) 

ROUT Sll)=RuUTS (2 )=ROGT S l 3) = T 

IF ( T.L T.O..GR. T.GT. 1.) GG TO 100 
GU TO 90 

TEMP=CGEFFS{ 3) **2-4. *CGEFF S ( 2 )*CGEFFS ( 4 ) 

IF l TEMP.LT.G. ) GG TO 10 
TMP=SWKT(TEMP) 

ROOT 3 ( 1)=RG0TS 12 ) = (- CGEFFSl 3 )+TMP)/ (2.*CGEFFS( 2J ) 
ROUT SI 3)-(-CGEFF SI 3 J-TMP I / (2. *COEFFS 12) ) 

GG TG 30 
CGNT INGE 

CALL CUBIC ICGEFFS, ROOTS) 


SELECT DESIRED ROOT 
C 

TT^REALIROGTSI i) ) 

IF ( AIMAGI RGGTSl 1) ) .EU.O. . ANO.TT .GE .0. . ANJi. TT. LE. 1. ) GU TG 60 
S»0 CONTINUE 

TT=REALlRG0TS(2) ) 

IF (AIMAGI ROGTSl 2) ) .Ew.O. .AND. TT .GE .0. . AND.. T T. LE . 1. ) GG TG BO 
TT=REAL(R00TS13) ) 

IF ( AIMAGI ROOT Si 3 )). EQ. 0. .AND .TT .GE .0. .AND.T T. LE. 1. ) GG TO 40 

GU TO LUO 

60 IF (AIM AG (ROOTS! 2) ) .NE. 0. ) GO TO 90 
TT=REALi ROOT S( 2) ) 

IF ( TT.LT.0..0R.TT.GT.1.) GO TO BO 
PRINT 70, T» TT , TL»CUtFFS»ROUT S 

70 FORMAT (//40H TWO ACCEPTABLE ROOTS FOUNG» ROUTS ARE *2E17.B/34H E 
ESTIMATE FUR SELECTION OF ROOTS ,E17.3/15H COEFFICIENTS =,4E17.8/8 
2H ROOTS -=|6E1 7.8// ) 

IF (ABS(TL-TT).LT.AfaS(TL-T)) T=TT 
60 IF l A I MAG( RGGTS13) ).NE.O.) GC TO 90 

TT=REAL(RG0TSl3i ) 

IF ( TT.LT.0..GR.TT.GT.1.) GO TO 90 
PRINT 70, T, TT»TL, CO EFFS, ROOTS 
IF ( ABSITL-T T) .LT. AfaSt TL-T) ) T=TT 
90 CONTINUE 

K U 8S GL = 1 
RETURN 

100 CONTINUE 
K.U6SUL=2 
RETURN 
END 
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Subroutine CUBIC 


Subroutine CUBIC uses the direct solution for the roots of the cubic equation: 

AX'* + BX^ + CX + D = 0. The flow chart and the FORTRAN statements for this subroutine 
are as follows: 


C 


Subroutine 

CUBIC 


5 

7 

Solve directly 
for the roots 
of a cubic 
equation 




C 


) 


SUBROUTINE CUBIC IA,X) 

C 

c 

C PROGRAMMER HSING-SHIH, CHANG 

C 

C 

COMPLEX X 
DIMENSION X(3) 

DIMENSION A(4I,XR13J, XIOJ.AUOI 
C 

C SOLVE THE CUBIC EQUATION F(X)=0 

C 

C F(X)=A(n*X**3+A(2)*X**2+A(3)*X+A(4) 

C 

C QUADRA T DON FIX)=0 

C 

C F ( X) =Aq( 1J*X**2*AQ( 2J*X+AQ(3) 

C 

C F(X)=X**2+B2*X+B3 

I PATH— 2 
EX- I./ 3. 

IF (AC4J) 20,10,20 
10 XR 1 1 ) =0 . 

GO TO 150 
20 A2=A(1J*A(11 

Q=t27.*A2*A{4)-9.*A( 1 ) * A ( 2 ) * A { 3 ) +2 . *A( 2 )**3) / ( 54.*A2*A(1) ) 
IF IQ) 40,30, 30 
30 2=0. 

GO TO 140 

<►0 u=-u 

I PAT H= 1 

30 P=(3.*AI1I *A ( 3 )- A( 2) #A( 2) )/ { 9.*A2 ) 

ARG=P*P*P-*-Q*Q 
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IF ( ARG1 60,70,80 

60 Z=-2.*SURT (-P)*CUS(ATAN( SORT (-ARG1 /Q 1 /3. ) 

GO TO 120 
70 Z=-2.*U**EX 

GO TO 120 

80 SARG=SURT( ARG) 

IF (PI 90 , 100,110 

90 Z=-(U«-SARG1**EX-(Q-SARG1**EX 

GO TO 120 

100 Z=-(2. *Q)**EX 

GU TO 120 

110 Z=(SARG-Q)**EX-(SARG*Q)**EX 
120 GO TO (130,140), IPATH 
130 L-—L 

140 XR(11=(3.*A( 1) *Z-A(2l )/ (3,*A( 1) ) 

130 AO ( 1 ) = A ( 1 1 

A0(21=A(21 *XR( il*A( 1) 

A0(31 = A(3H-XR(ll*A0t21 
C 

62=AQ(21/AQ( 11 
B3=AQ(3)/AQ( 1) 

Xl=-62/2. 

OISC=Xl*Xl— 83 
IF (OISC.LT.O.O) 160,170 
160 X2=SQRT (-DISC) 

XR(21=X1 
XR ( 3 1 =X1 
XX (2)=X2 
GO TO 200 

i/0 IF (DISC. EU. 0.01 180,190 

160 X2=0 . 

XR(2 )=X1+X2 
XR ( 3 1 = X1-X2 
Xl(2)= 0. 

GO 10 200 

190 X2=SURT(U1 30 

XR(21=X1+X2 
XR(3)=X1-X2 
XI (21=0. 

200 XX(il=0. 

Xl(3)=-XI(2J 
x(u=cmplx(aR( li ,xk m 
X(2) =CMPLX(XR(2) ,XI(21 ) 

X ( 31 =CMPCX ( XR ( 3 1 , X I ( 3 1 1 

RETURN 

ENO 



Subroutine PLTROT 


Subroutine PLTROT rotates and translates points defining a cross section, generates 
plot instructions, and prints the points. The flow chart and the FORTRAN statements for 
this subroutine are as follows: 


C 


Subroutine PLTROT 


) 


Print 
original 
line of 
points 


Rotate and 
project 
points 



r 


_3Z_ 


Print 
/ rotated line 

[ ? f 

V points 


W 

^ RETURN ^ 


SUBROUTINE PLTROT (NPT, A, ALINE, HMIN, VMlNtSCALE » IPRUNtXP, NAME) 

c 

C ROTATES A SET OF 30 POINTS INTO A SPECIFIED PLANE 

C AND GENcRATES A CALCCMP PLOT TAPE 

C 

DIMENSION A(3»3) »ALINE(52»3) »RLINE(54»3) 

DIMENSION XP 1 3 ) f NAME ( 2 1 
DATA EPS/. 00000001/ 

c 

N=1 

10 N=N+1 

Tl=A8S(ALINE(N-l,l)-ALlNE(N, 11) 

T2=A6S(ALIN£(N-1 t 2I-AL I NE ( N » 2 ) ) 

T3=A6S (ALINE (N— 1 » 3 1 -AL I NE ( N* 3) ) 

IF (. NOT. (Tl. LE.EPS.AN u.T 2.LE.EPS. ANO.T3.LE.EPS)! GO TO 30 
DQ 20 NK=N» NPT 
DO 20 N3=l »3 

20 ALINE ( NK-I »N3J=ALINE(NK»N3I 

NPT=NPT-i 
N=N— 1 
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30 IF (N. WE. NPT) 00 TO 10 

IF INPT.LE.il RETURN 
GO TO (70,40,70,40) , IPR1N 
40 CONTINUE 

00 50 N=1 » NP T 
DO 50 N3=l,3 

30 RLiNE(N»N3 )= ALInc( N, N3) /SCALE 

WRITE (6,60) ( I ALINE ( I » J) » J= l»3) » (RLINEII . J) » J=l* 3) • 1 = 1, NPT) 
oG FORMAT ( //44X, 3/HCOORul NATES OF POINTS OF INTERSEC T I ON/2 7X , 8H0R IGI 
INAL52X,6HSCALEG//7X, 1HX 19X , 1HY 19X, 1 H*. 1 9X, 1 HX 19X, IHY 19X, lHZ/( 6E20.8 

2) ) 

70 CONTINUE 

DO 100 N= 1 »NPT 
RLINEINt 11*0. 

RL INE ( N» 2 ) =0 . 

RL1NE(N,3)=U. 

DO 90 1=1,3 
DO 80 J= l » 3 

60 RLI.NE(N,I)=RLINE(N,i)-*-A( I , J i *AL INE { N, J ) 

90 KL INE(N, I )=RLINE(N, l )-XP( l I 

100 CONTINUE 

RLINE(NPT*1»2) =HM I N 
RLINE(NPT*1,3)=VMIN 
RLINE(NPT*2»2)=SCALE 
KLINE I NPT *2,3) =SCALE 

CALL LINE (RLINEI 1,2) ,RLINE( 1, 3) ,NPT, 1 ,0,0,0) 

GO TO (150,150,110,110), IPRIN 
110 CONTINUE 

DO 130 N=1,NPT 
DO 120 N3= 1,3 

120 AL 1NE( N , N3 )=RL I NE ( N, N3 ) /SCALE 

130 CONTINUE 

WRITE (6,140) ( (RL INEII , J) , J=l,3) , ( ALINEH ,J) , J=l, 3) , 1 = 1, NPT) 

140 FORMAT ( // 33 X ,59HC00RDI NATES OF POINTS OF INTERSECTION ROTATED INT 
10 Y L PLANE/2 7X , 6HQR I G INAL52X , 6HSCAL ED/ / 7X, 1HX19X, 1HY19X, 1HZ19X, IHX 
219X, IHY19X, l HZ/ (oE 20. 8) ) 

150 CONTINUE 
RETURN 

C 1 

C END OF PLTRUT 

C 

END 
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Function IPCOF 


Function IPCOF finds the coefficients of a plane from three given points, 
chart and the FORTRAN statements for this function are as follows: 


C 


Function IPCOF 


) 


± 

Compute the 
equation of a 
plane from 
3 given points 


r^L 


5 


FUNCUCN IPCOF I A, B,C,0,81,P2,P3 J 

DIMENSION PllJi ,P2(3)»P3(3)»AMAT(3»3) ,BMAT(3J , IPIVGT13) 
C 

C FINDS THE COEFFICIENTS OF A PLANE IN THE 

C FORM AX+bY+CZ=D 

C 

DATA EPS/1. E-b/ 

C 

0=1. 

DO SO 1=1,4 
DO 10 N 3=1,0 
AM AT 1 1 ,N3)=P1 (Nil 
AM AT 1 2 ,N3 J =P2( N3 ) 

AMAT ( 3 ,N3) =P3(N3 ) 

10 BMAT I N31 = i • 

11=1-1 

IF III. EO. 01 oO TO 30 
DO 20 N3=l , 3 
bMAT (N3J=-AMAT IN 3, III 
AMAT ( N 3 , 1 1 i =— 1 . 

20 CONTINUE 

30 CONTINUE 

CALL SIMEQ I AMAT ,3 , OMAT , 1 ,U£T , IP I VOT , 3 , 1SCALE J 
IF l AbSlOET) .Lt.EPSl GO TO 50 
IF I l.EQ.l 1 GO TO 40 
D= BMAT 1 1 l ) 

8MAT ( 1 1 I = 1 . 

40 A=BMAT 1 1 1 

B= BMAT 121 
C=BMAT(3J 
I PC0F= 1 
RETURN 

50 CONTINUE 

I PC0F=2 
RETURN 
C 


C 

END 


The flow 
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Langley Library Subroutine SIMEQ 


Language : FORTRAN 

Purpose: SIMEQ solves the matrix equation AX = B where A is a square coefficient matrix and B is a 

matrix of constant vectors. The solution to a set of simultaneous equations and the determinant may be 
obtained. If the user wants the determinant only, use DETEV for savings in time and storage. 

Use : CALL SIMEQ (A, N, B, M, DETERM, IPIVOT, NMAX, ISCALE) 

A A two-dimensional array of the coefficients. 

N The order of A; 1 i N ^ NMAX. 

B A two-dimensional array of the constant vectors B. On return to calling program, X is 

stored in B. 

M The number of column vectors in B. 

DETERM Gives the value of the determinant by the following formula: 

DET(A) = (io 100 ) ISCALE (DETERM) 

IPIVOT A one -dimensional array of temporary storage used by the routine. 

NMAX The maximum order of A as stated in dimension statement of calling program. 

ISCALE A scale factor computed by subroutine to keep results of computation within the floating- 
point word size of the computer. 

Restrictions: Arrays A, B, and IPIVOT are dimensioned with variable dimensions in the subroutine. The 

maximum - size of these arrays must be specified in a DIMENSION statement of the calling program as: 

A (NMAX, NMAX), B (NMAX, M), IPIVOT (NMAX). The original matrices, A and B, are destroyed. 

They must be saved by the user if there is further need for them. The determinant is set to zero for 
a singular matrix. 

Method: Jordan’s method is used through a succession of elementary transformations: Z n , ^n-l> ■ • *» *1’ 

If these transformations are applied to a matrix B of constant vectors, the result is X where AX = B. 
Each transformation is selected so that the largest element is used in the pivotal position. 

Accuracy: Total pivotal strategy is used to minimize the rounding errors; however, the accuracy of the 
final results depends upon how well -conditioned the original matrix is. 

Reference: (a) Fox, L.: An Introduction to Numerical Linear Algebra. Oxford Univ. Press, c.1965. 

Storage : 432 q locations. 

Subroutine date: August 1, 1968. 
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The FORTRAN statements for this subroutine are as follows: 

SUBROUTINE S IMEQ ( A, N , B » M , OE T ERM , I P I VQT ,NMAX , I SC AL E J 
C SOLUTION OF SIMULTANEOUS LINEAR EQUATIONS 

C *** DOCUMENT DATE 06-01-68 SUBROUTINE REVISED 08-01-68 ********* 

C 

DIMENSION IP l VOT (N J , AIN MAX , N I ,B(NMAX,M) 

EQUIVALENCE ( IRQ*, JRGWJ , ( I COLUM , JCOLUM 1 , ( AMAX , T , SW AP i 
C 

C INITIALIZATION 

c 

10 ISCAlE=0 

Rl=lO. 0**100 
R2=l.0/Rl 
DET£RM=1 .0 
DO 20 J= l » N 
20 IPIVOT ( J) = 0 

DO 380 1 = 1, N 
C 

C SEARCH FOR PIVOT ELEMENT 

C 

AMAX=0.0 
DO 70 J=1,N 

IF ( IPIVOT 1 J )-I ) 30,70,30 
30 DO 60 K=1,N 

IF ( IP IVOT t K I — II 40,60,300 
40 IF I ABSl AM AX )— ABSI A ( J »K ) I ) 30,60,60 
30 1R0W=J 

ICOLUM=K. 

AMAX=A( J ,K> 
oO CONTINUE 

70 CONTINUE 

IF (AMAX) 90,80, 90 
80 DtTERM=0.0 

1SCALE=0 
GO TO 390 

90 IPIVOT ( ICOLUM)= I PI VOT ( I CCLUM ) +1 

C 

C INTERCHANGE ROWS TO PUT PIVOT ELEMENT ON DIAGONAL 

C 

IF ( I ROW— I COLUM) 100,140,100 
100 DETERM=— DETERM 
00 110 L= 1 , N 
SWAP=A 1 1 ROW, L ) 

A ( 1K0R ,L ) = A 1 ICO(.UM,L ) 

110 A( 1C0LUM.L )=SWAP 

IF (M) 140,140,120 
120 DO 130 L=l,M 

SWAP=B( IROW»L) 
b( 1R0W »L)=6( ICOLUM,L) 

130 81 ICOLUM,L)=SWAP 

140 PI VGT=A ( 1 COLUM, I COLUM) 

IF (PIVOT) 150,80,130 
C 

C SCALE THE DETERMINANT 

C 

150 PIVOT I=P I VOT 

IF (ABS(DETERM)-Rl ) 180,160,160 
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160 

170 

180 

150 


200 

210 

220 


<:30 


240 

230 


260 

270 

C 

C 

C 


280 

290 

300 

310 

C 

c 

c 

320 

330 


340 

350 

360 

370 

380 

390 


DET£RM=UETERM/kl 
ISCALE = ISCAlE«- 1 
IF ( ABSIOETERMJ-RU 
0c TERM=DtT ERM/R1 
1 SCALE- I SC AL E+ l 
GO TO 210 

IF IA8SIDETERM)-R2J 
UE T£RM=OET ERM*R1 
l SCA lE= ISC AL E- 1 
IF IA8SIDETERMJ-R2) 
DETERM=uETERM*Ki 
I SCalE=ISCALE- 1 
IF lAdSIPIVOTI) — Rl) 
PIVOTl=PlVUT I/Rl 
ISCAi_E=lSCALE«-l 
IF IA6SIPI VUTIJ-RU 
PIVOT I=P IVGT l/Rl 
ISCALE=ISCALE+1 
GO TO 270 

IF (A8SIPIVUTD-R2) 
PIVOT I=PIVOT l*kl 
ISCAL6=ISCAL£-1 
IF (A6SIPIV0T I )-K2J 
PIVOT l=PIVOTI*Rl 


210,170,170 

190. 150.210 

200 . 200.210 
240,220,220 
270,230,230 

250 .250.270 

260.260.270 


1SCALE=ISCALE-1 
UETERM=0ETERM*P1 VOTI 


DIVIDE PIVOT RUrt 8 Y PIVOT ELEMENT 
DO 250 L=1,N 

IF I 1 PIVOT CL 1 — 1 J 280,290,390 
A I 10 UL UM » L ) — A ( ICCLUM.L) /PIVOT 
CONT INOE 

IF (M) 320,320,300 
DO 310 L=1,M 

b I ICOLUM,L)=Bl ICOLUM.L) /P1VGT 


REDUCE NON-PIVOT RUrtS 

DO 380 L 1= 1 , N 
IF ILl-lCOLUMl 330,380,330 
T=AILl,I CULUM J 
DU 350 L=1 ,N 

IF (lPIVUTCL)-n 340,350,390 
A(L1,L)-A(L1,L)-At ICOLUM,LJ*T 
CONTINUE 

IF (M) 380,380,360 
DU 370 L=1,M 

8(LI,U=BIL1,U-B( I COLUM, L ) *T 

CONTINUE 

RETURN 

END 



PROGRAM USE 


PROGRAM IDENTIFICATION 

This program is for fitting smooth surfaces to the component parts of an aircraft 
configuration using a three-dimensional modeling technique called Coon's patches (ref. 1). 
It is identified as program D3400. 

PROGRAM SETUP FOR A COMPILE AND EXECUTE 

This section describes the input data requirements, limitations, and the punched 
card formats. The program will end normally if there are no input cards at the beginning 
of a READ sequence. 

The input data cards are assembled with the program decks in the order illustrated 
in sketch (d). 
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DESCRIPTION OF INPUT DATA CARDS 


Configuration 

The form for the airplane configuration input has become known throughout the air- 
craft industry as the Harris Wave Drag geometry input and is identical to that described 
in reference 3. 

Since the airplane has to be symmetrical about the XZ-plane, only half of the air- 
plane heed be described to the computer. The convention used in presenting the input data 
is that the half of the airplane on the positive Y-side of the XZ-plane is presented. The 
program then uses this information to construct the complete airplane if required. The 
number of input cards depends on the number of components used to describe the configu- 
ration and the amount of detail used to describe each component. It is not possible to 
change only part of a configuration in a succeeding case. The complete configuration 
must be input each time. The method of input is by FORTRAN "READ" statements. 

Card 1 - Identification.- Card 1 contains any desired identifying information in 
columns 1 to 80. 

Card 2 - Control integers .- Card 2 contains 24 integers, each punched right- 
justified in a 3-column field. Columns 73 to 80 may be used in any desired manner. An 
identification of the card columns, the name used by the source program, and a descrip- 
tion of each integer are given in the following table: 


Columns 

FORTRAN 

name 


Description 

01 to 03 

JO 

If 

JO = 0, no reference area 



If 

JO = 1, reference area to be read 

04 to 06 

J1 

If 

Jl = 0, no wing data 



If 

Jl = 1, cambered wing data to be read 



If 

Jl = -1, uncambered wing data to be read 

07 to 09 

J2 

If 

J2 = 0, no fuselage data 



If 

J2 = 1, data for arbitrarily shaped fuselage 




to be read 



If 

J2 = -1, data for circular fuselage to be read 


(with J6 = 0, fuselage will be 
cambered; with J6 = -1, fuselage 
will be symmetrical with XY-plane; 
with J6 = 1, entire configuration 
will be symmetrical with XY-plane) 
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Columns 

FORTRAN 

name 



Description 

10 to 12 

J3 

If 

J3 = 0, 

no pod data 



< If 

J3 = 1, 

pod data to be read 

13 to 15 

J4 

If 

o' 

II 

no fin data 



If 

J4 = 1, 

fin data to be read 

16 to 18 

J5 

If 

J5 = 0, 

no canard data 



If 

J5 = 1, 

canard data to be read 


19 to 21 

J6 

Simplification code: 

If J6 = 0, indicates a cambered circular or 
arbitrary fuselage if J2 ^ 0 
If J6 = 1, complete configuration is symmet- 
rical with respect to XY -plane, 
which implies uncambered circular 
fuselage if there is a fuselage 
If J6 = -1, indicates uncambered circular 
fuselage with J2 ^ 0 

22 to 24 

NWAF 

Number of airfoil sections used to describe the 
wing; 4 = NWAF = 20 

25 to 27 

NWAFOR 

Number of ordinates used to define each wing 
airfoil section; 4 s NWAFOR s 30 

28 to 30 

NFUS 

Number of fuselage segments; 1 £ NFUS ^ 4 

31 to 33 

NRADX(l) 

Number of points used to represent half -section 
of first fuselage segment; if fuselage is circu- 
lar, the program computes indicated number of 
y- and z-ordinates; 4 g NRADX(l) S 30 

34 to 36 

NFORX(l) 

Number of stations for first fuselage segment; 
4 g NFORX(l) ^ 30 

37 to 39 
40 to 42 

NRADX(2) 

NFORX(2) 

Same as NRADX(l) and NFORX(l), but for second 
fuselage segment 

43 to 45 
46 to 48 

NRADX(3) 

NFORX(3) 

Same as NRADX(l) and NFORX(l), but for third 
fuselage segment 

49 to 51 
52 to 54 

NRADX(4) 

NFORX(4) 

Same as NRADX(l) and NFORX(l), but for fourth 
fuselage segment 
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Columns 

FORTRAN 

name 

Description 

55 to 57 

NP 

Number of pods described; NP ^ 9 

58 to 60 

NPODOR 

Number of stations at which pod radii are to be 
specified; 4 s NPODOR S 30 

61 to 63 

NF 

Number of fins (vertical tails) described; 
NF ^ 6 

64 to 66 

NFINOR 

Number of ordinates used to define each fin air- 
foil section; 4 S NFINOR i 10 

67 to 69 

NCAN 

Number of canards (horizontal tails) described 
NCAN i 2 

70 to 72 

NCANOR 

Number of ordinates used to define each canard 
airfoil section; 4 ^ NCANOR ^ 10; if NCANOR 
is given a negative sign, the program will expect 
to read lower ordinates also; otherwise, airfoil 
is assumed to be symmetrical 


Car ds 3, 4, ... - remaining data input cards .- The remaining data input cards 
contain a detailed description of each component of the airplane. Each card contains up 
to 10 values, each value punched in a 7-column field with a decimal and may be identified 
in columns 73 to 80. The cards are arranged in the following order: reference area, 
wing data cards, fuselage data cards, pod (or nacelle) data cards, fin (vertical tail) data 
cards, and canard (or horizontal tail) data cards. 

Reference area card: The reference area value is punched in columns 1 to 7 and 

may be identified as RE FA in columns 73 to 80. This value is not used by the program 
but may be present in an already existing input deck. 

Wing data cards: The first wing data card (or cards) contains the locations in per- 
cent chord at which the ordinates of all the wing airfoils are to be specified. There will 
be exactly NWAFOR locations in percent chord given. Each card may be identified in 
columns 73 to 80 by the symbol XAFj where j denotes the number of the last location in 
percent chord given on that card. For example, if NWAFOR = 16, there are 16 ordinates 
to be specified for every airfoil, and two data cards will be required. The first XAF card 
is identified as XAF 10 and the second as XAF 16. 

The next wing data cards (there will be NWAF cards) each contain four numbers 
which give the origin and chord length of each of the wing airfoils that is to be specified. 
The cards representing the most inboard airfoil are given first, followed by the cards for 
successive airfoils. The information is arranged on each card as follows: 
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Columns 


Description 


1 to 7 
8 to 14 
15 to 21 
22 to 28 
73 to 80 


x-ordinate of airfoil leading edge 

y -ordinate of airfoil leading edge 

z -ordinate of airfoil leading edge 

Airfoil streamwise chord length 

Card identification, WAFORGj where j denotes the 
particular airfoil; for example, WAFORG1 denotes 
first (most inboard) airfoil 


If a cambered wing has been specified, the next set of wing data cards is the mean 
camber line (TZORD) cards. The first card contains up to 10 Az values, referenced 
to the z-ordinate of the airfoil leading edge, at each of the specified percents of chord 
for the first airfoil. If more than 10 values are to be specified for each airfoil (there 
will be NWAFOR values), the remaining values are continued on successive cards. The 
remaining airfoils are described in the same manner, data for each airfoil starting on a 
new card, and the cards arranged in the order which begins with the most inboard airfoil 
and proceeds to the outboard. Each card may be identified in columns 73 to 80 as 
TZORDj, where j denotes the particular airfoil. 

Next are the wing airfoil ordinate (WAFORD) cards. The first card contains up to 
10 half -thickness ordinates of the first airfoil expressed as percent chord. If more than 
10 ordinates are to be specified for each airfoil (there will be NWAFOR values), the 
remaining ordinates are continued on successive cards. The remaining airfoils are each 
described in the same manner, and the cards are arranged in the order which begins with 
the most inboard airfoil and proceeds to the outboard. Each card may be identified in 
columns 73 to 80 as WAFORDj, where j denotes the particular airfoil. 

Fuselage data cards: The first card (or cards) specifies the x-values of the fuse- 
lage stations of the first segment. There will be NFORX(l) values and the cards may be 
identified in columns 73 to 80 by the symbol XFUSj where j denotes the number of the 
last fuselage station given on that card. 

If the fuselage is circular and cambered, the next set of cards specifies the 
z-locations of the center of the circular sections. There will be NFORX(l) values and 
the cards may be identified in columns 73 to 80 by the symbol ZFUSj where j denotes 
the number of the last fuselage station given on that card. 

If the fuselage is circular, the next card (or cards) gives the fuselage cross- 
sectional areas, and may be identified in columns 73 to 80 by the symbol FUSARDj where 
j denotes the number of the last fuselage station given on that card. If the fuselage is of 
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arbitrary shape, the y-ordinates for a half -section are given (NRADX(l) values) and 
identified in columns 73 to 80 as Yi where i is the station number. Following these 
are the corresponding z -ordinates (NHADX(l) values) for the half -section identified in 
columns 73 to 80 as Zi where i is the station number. Each station will have a set of 
Y and Z cards, and the convention of ordering the ordinates from bottom to top is 
observed. 

For each fuselage segment a new set of cards as described must be provided. The 
segment descriptions should be given in order of increasing values of x. 

Pod data cards: The first pod or nacelle data card specifies the location of the 
origin of the first pod. The information is arranged on the card as follows: 

Columns Description 

1 to 7 x-ordinate of origin of first pod 

8 to 14 y-ordinate of origin of first pod 

15 to 21 z-ordinate of origin of first pod 

73 to 80 Card identification, PODORGj where j denotes pod number 

The next pod input data card (or cards) contains the x-ordinate s, referenced to the 
pod origin, at which the pod radii (there will be NPODOR of them) are to be specified. 

The first x-value must be zero, and the last x-value is the length of the pod. These cards 
may be identified in columns 73 to 80 by the symbol XPODj where j denotes the pod 
number. For example, XPOD1 represents the first pod. 

The next pod input data cards give the pod radii corresponding to the pod stations 
that have been specified. These cards may be identified in columns 73 to 80 as PODRj 
where j denotes the pod number. 

For each additional pod, new PODORG, XPOD, and PODR cards must be provided. 
Only single pods are described but the program assumes that if the y-ordinate is not zero 
an exact duplicate is located symmetrically with respect to the XZ-plane; a y-ordinate of 
zero implies a single pod. 

Fin data cards: Exactly three data input cards are used to describe a fin. The 
information presented on the first fin data input card is as follows: 

Description 

x-ordinate of lower airfoil leading edge 
y-ordinate of lower airfoil leading edge 
z-ordinate of lower airfoil leading edge 
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Columns 
1 to 7 
8 to 14 
15 to 21 



Columns 

Description 

22 to 28 

Chord length of lower airfoil 

29 to 35 

x-ordinate of upper airfoil leading edge 

36 to 42 

y-ordinate of upper airfoil leading edge 

43 to 49 

z -ordinate of upper airfoil leading edge 

50 to 56 

Chord length of upper airfoil 

73 to 80 

Card identification, FINORGj where j denotes fin number 


The second fin data input card contains up to 10 locations in percent chord (exactly 
NFINOR of them) at which the fin airfoil ordinates are to be specified. The card may be 
identified in columns 73 to 80 as XFINj where j denotes the fin number. 

The third fin data input card contains the fin airfoil half -thickness ordinates 
expressed in percent chord. Since the fin airfoil must be symmetrical, only the ordinates 
on the positive y-side of the fin chord plane are specified. The card identification, 
FINORDj, may be given in columns 73 to 80, where j denotes the fin number. 

For each fin, new FINORG, XFIN, and FINORD cards must be provided. 

Only single fins are described but the program assumes that if the y-ordinate is not 
zero an exact duplicate is located symmetrically with respect to the XZ-plane; a 
y-ordinate of zero implies a single fin. 

Canard data cards: If the canard (or horizontal tail) airfoil is symmetrical, 
exactly three cards are used to describe a canard, and the input is given in the same 
manner as for the fin. If, however, the canard airfoil is not symmetrical (indicated by a 
negative value of NCANOR), a fourth canard data input card will be required to give the 
lower ordinates. The information presented on the first canard data input card is as 


follows: 


Columns 

Description 

1 to 7 

x-ordinate of inboard airfoil leading edge 

8 to 14 

y-ordinate of inboard airfoil leading edge 

15 to 21 

z -ordinate of inboard airfoil leading edge 

22 to 28 

Chord length of inboard airfoil 

29 to 35 

x-ordinate of outboard airfoil leading edge 

36 to 42 

y-ordinate of outboard airfoil leading edge 

43 to 49 

z -ordinate of outboard airfoil leading edge 
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Columns 


Description 


50 to 56 Chord length of outboard airfoil 

73 to 80 Card identification, CANORGj where j denotes the canard number 

The second canard data input card contains up to 10 locations in percent chord 
(exactly NCANOR of them) at which the canard airfoil ordinates are to be specified. The 
card may be identified in columns 73 to 80 as XCANj where j denotes the canard 
number. 

The third canard data input card contains the upper half -thickness ordinates, 
expressed in percent chord, of the canard airfoil. This card may be identified in 
columns 73 to 80 as CANORDj where j denotes the canard number. If the canard air- 
foil is not symmetrical, the lower ordinates are presented on a second CANORD card. 
The program expects both upper and lower ordinates to be punched as positive values in 
percent chord. 

For another canard, new CANORG, XCAN, and CANORD cards must be provided. 

Alternate Surface -Description Input 

The surface -description method used in this report is not restricted to any partic- 
ular shape. An airplane configuration was chosen for the current application because of 
its complexity and an anticipated need. The program input is best thought of as a set of 
points describing a surface or surfaces. 

To use the program for a different input form, the user may substitute another 
overlay (1,0) for preprocessing input data in the proper form to be used by the program. 
If desired to print this input data, provisions to do so should be provided for in over- 
lay (1,0). 

Overlay (1,0) must: 

1. Store in labeled COMMON 

COMMON /PATPLT /XMIN, XMAX, YMIN, YMAX, ZMIN, ZMAX, NOBJ 

The minimums and maximums define a box in which the shape to be plotted ortho- 
graphically will fit. These values are not used for the cross-section plots, but space 
must be allowed so that NOBJ will occupy the correct position. 

NOBJ is the number of objects (or components) each made up of a number of sur- 
faces, all of which could form a body. 

2. Write binary tape 10 in the following format. FORTRAN names from the given pro- 

gram are used for illustration. 
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Record 


Name 


Purpose 



ABC(8) 

Identification 

(For NOBJj) 


NSURF 

Number of surfaces in object 

Ml 

Not used 

M2 (2) 

Name of object (for printing only) 

M3 

Not used 

M4 

Not used 

(For NSURF -jJ 


NCOL 

Number of columns of grid points, ^31 

NROW 

Number of rows of grid points, ^31 

N3 

Not used 

N4 

Not used 

N5 

Not used 

((ALENE(N,N3), N=l, NCOL), 
N3=l,3) 

NROW lines of NCOL points (x,y,z) 

((ALINE(N,N3), N=l, NROW), 

NCOL lines of NROW points (x,y,z) 


N3=l,3) 

(Repeat records 3, 4, and 5 for the number of surfaces given in record 2) 


(Repeat record 2 for each NOBJ as given in labeled COMMON PATPLT 
and repeat records 3, 4, and 5 as required.) 

(Although the dummy variables are not used at the present time, they 
must be written on tape.) 


Great care must be taken to describe the grid of input points in exactly the manner 
specified so that in a collection of surfaces the outward normal vectors will be consistent. 
For the current application, the rectangular grid of values (not necessarily rectangular in 
shape) does not require any specific corner as the starting place. However, points must 
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be given as lines in a rowwise direction then as lines in a columnwise direction as 
illustrated in sketch (e) 



3. If it is desired to increase the number of points in the rows and/or columns and it is 
not feasible to break a surface down into more than one surface, increase all dimen- 
sions of 31 to the maximum of NCOL or NROW in: 

(a) Overlay (1,0), program START 

(b) Overlay (2,0), subroutine PACH. Also in subroutine PACH, change the dimen- 

sions of the SLOPE array to SLOPE(NCOL max ,NROW max ,3) and change 
MAXN in the data statement to the largest of NCOL max or NROW max - 


Option Card 

The option card indicates to the program the next kind of input to be read. 


Column 

4 


FORTRAN 

name 


ITYPE 


Description 

If ITYPE = 1, Read orthographic-projection plot 
card 

If ITYPE = 2, Read cross-section plot card 
If ITYPE = 3, Read another set of geometry cards 


Plot Cards 

A single card contains all the necessary information for one plot. The available 
options and the necessary input for each are described in the following sections. 


Orthographic projections (ITYPE = 1) .- For one orthographic projection, the card 
should be set up as follows: 


Columns 


FORTRAN 

name 


Description 


1 

HORZ 

"X", "Y", or "Z" for horizontal axis 

3 

VERT 

"X", "Y", or "Z" for vertical axis 

5 to 7 

TEST1 

Word "OUT" for deletion of hidden lines; otherwise, 
leave blank 

8 to 12 

PHI 

Roll angle, degrees 

13 to 17 

THETA 

Pitch angle, degrees ' ' 

18 to 22 

PSI 

Yaw angle, degrees 

48 to 52 

PLOTSZ 

Plot frame size (scale factor is computed by using 
PLOTSZ and maximum dimension of configuration, 
i.e., the maximum dimension, usually the body 
length, will be scaled to exactly PLOTSZ) 

53 to 55 

TYPE 

Word "ORT" 

56 to 58 

NOU 

Number of lines, originating on the w = 0.0 boundary, 
computed within each patch for enriching the sur- 
face grid in the w-direction, NOU ^ 50 

59 to 61 

NOW 

Number of lines, originating on the u = 0.0 boundary, 
computed within each patch for enriching the sur- 
face grid in the u-direction, NOW s 50 

64 

ISIDE 

If ISIDE = 0 or 1, plot described body 

If ISIDE = 2, plot described body and mirror image 

72 

KODE 

If KODE = 0, continue to read plot cards 
If KODE ^ 0, return and read option card after this 


plot 


Plan, front, and side views (stacked) (ITYPE = 1) .- For plan, front, and side views 
stacked one above the other in a pleasing -to-the-eye arrangement, the card should be set 
up as follows: 


Columns 


FORTRAN 

name 


Description 


5 to 7 TEST1 Word "OUT” for deletion of hidden lines; other- 

wise leave blank 
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Columns 


FORTRAN 

name 


Description 


8 to 12 
13 to 17 
18 to 22 
48 to 52 


53 to 55 
56 to 58 


PHI 

THETA 

PSI 

PLOTSZ 


TYPE 

NOU 


59 to 61 NOW 


64 ISIDE 


72 


KODE 


y -origin on paper of plan view, inches 

y -origin on paper of side view, inches 

y -origin on paper of front view, inches 

Plot frame size (scale factor is computed 
using PLOTSZ and maximum dimension of 
configuration, i.e., the maximum dimension, 
usually the body length, will be scaled to 
exactly PLOTSZ) 

Word "VU3" 

Number of lines, originating on the w = 0.0 
boundary, computed within each patch for 
enriching the surface grid in the w -direction, 
NOU ^ 50 

Number of lines, originating on the u = 0.0 
boundary, computed within each patch for 
enriching the surface grid in the u -direction, 
NOW s 50 

If ISIDE = 0 or 1, plot described body 
If ISIDE = 2, plot described body and mirror 
image 

If KODE = 0, continue to read plot cards 
If KODE ± 0, return and read option card after 
this plot 


Cross-section plots where the method of input for defining the plane is by three 
input points (ITYPE = 2) .- For cross-section plots where the method of input for defining 
the plane is by three input points (ITYPE = 2), the input card should be set up as follows: 


Columns 

FORTRAN 

name 


Description 

1 to 6 

PPL1(1) 

xj, also used as 

x 0 

7 to 12 

PPL1(2) 

yj, also used as 

yo 

13 to 18 

PPL1(3) 

zj, also used as 

z 0 
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Columns 

FORTRAN 

name 

Description 

19 to 24 

PPL2(1) 

x 2 

25 to 30 

PPL2(2) 

?2 

31 to 36 

PPL2(3) 

z 2 

37 to 42 

PPL3(1) 

x 3 

43 to 48 

PPL3(2) 

^3 

49 to 54 

PPL3(3) 

z 3 

55 to 60 

PLOTS Z 

Scale factor, inches/unit 

61 to 63 

HPAGE 

Horizontal paper origin, inches (a value of 0.0 



will overlay the plots) 

64 to 66 

VPAGE 

Vertical paper origin will be half of VPAGE, 



inches 

67 to 69 

INP 

Word "PNT" 

70 to 72 

NOU 

Number of u-values between 0.0 and 1.0, together 



with the u = 0.0 and u = 1.0 values, where 
interpolated values of w are to be computed, 
NOU ^ 50 

73 to 75 

NOW 

Number of w-values between 0.0 and 1.0, together 



with the w = 0.0 and w = 1.0 values, where 
interpolated values of u are to be computed, 
NOW i 50 

76 to 77 

*ICUT 

Number of additional cross-section plots desired 

*For ICUT /0, a second input card is required: 

Columns 

FORTRAN 

name 

Description 

1 to 6 

DX 

AX for Xj, x 2 , and Xg 

7 to 12 

DY 

AY for y v y 2 , and y g 

13 to 18 

DZ 

AZ for Zp z 2 , and Zg 

21 

IH 

If IH = 0, plots will be overlaid with HPAGE 


applied only after last plot 

If IH ^ 0, normal HPAGE spacing between each 
plot 
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Columns 


FORTRAN 

name 


Description 


78 


79 


80 


ISIDE 

If 

ISIDE = 0 or 1, examine given body 


If 

ISIDE = 2, examine given body and its mirror 
image 

IPRIN 

If 

IPRIN = 0 or 1, do not print intersection 
points 


If 

IPRIN = 2, print original and original rotated 
points of intersection 


If 

IPRIN = 3, print scaled and scaled rotated 
points of intersection 


If 

IPRIN = 4, print points of intersection- 

original and original rotated, 
scaled and scaled rotated 

KODE 

If 

KODE = 0, continue to read plot cards 


If 

KODE ^ 0, return and read option card after 
this plot 


Cross-section plots where the method of input for defining the plane is by specifying 
a Mach number to define a Mach angle, an angle to define the orientation of the Mach 
angle, and the plane intercept on the X-axis (ITYPE = 2) .- For cross-section plots where 
the method of input for defining the plane is by specifying a Mach number to define a Mach 
angle, an angle to define the orientation of the Mach angle, and the plane intercept on the 
X-axis (ITYPE = 2), the input card should be set up as follows: 


Columns 

FORTRAN 

name 

Description 

1 to 6 

PPL1(1) 

X 0 

7 to 12 

PPL1(2) 

Y 0 

13 to 18 

PPL1(3) 

z o 

19 to 24 

PPL2(1) 

X-intercept 

25 to 30 

PPL2(2) 

Roll angle of plane, degrees 

31 to 36 

PPL2(3) 

Mach number, >1.0 

55 to 60 

PLOTS Z 

Scale factor, inches/unit 


61 to 63 

HPAGE 

Horizontal paper origin, inches (a value of 0.0 
will overlay the plots) 

64 to 66 

VPAGE 

Vertical paper origin will be half of VPAGE, 
inches 
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Columns 

67 to 69 
70 to 72 


FORTRAN 

name 

INP 

NOU 


73 to 75 NOW 


76 to 77 *ICUT 

78 ISIDE 


79 IPRIN 


80 KODE 


Description 


Word "ANG" 

Number of u-values between 0.0 and 1.0, together 
with the u = 0.0 and u = 1.0 values, where 
interpolated values of w are to be computed, 
NOU ^ 50 

Number of w-values between 0.0 and 1.0, together 
with the w = 0.0 and w = 1.0 values, where 
interpolated values of u are to be computed, 
NOW ^ 50 

Number of additional cross-section plots desired 

If ISIDE = 0 or 1, examine given body 

If ISIDE = 2, examine given body and its mirror 
image 

If IPRIN = 0 or 1, do not print intersection points 

If IPRIN = 2, print original and original rotated 
points of intersection 

, If IPRIN = 3, print scaled and scaled rotated 
points of intersection 

If IPRIN = 4, print points of intersection: 

original and original rotated, 
scaled and scaled rotated 

If KODE = 0, continue to read plot cards 

If KODE / 0, return and read option card after 
this plot 


*For ICUT # 0, a second input card is 


Columns 

FORTRAN 

name 


1 to 6 

DX 

A 

7 to 12 

DY 

A 

13 to 18 

DZ 

A 

21 

IH 

If 


If 


required: 

Description 

for X-intercept 
for roll angle 
for Mach number 

EH = 0, plots will be overlaid with HPAGE 
applied only after last plot 

EH 1 0, normal HPAGE spacing between each 
plot 
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DESCRIPTION OF PROGRAM OUTPUT 


The program output includes the input data printout, the printout of points of inter- 
section for cross-section plots, and a plot vector file to be postprocessed for off-line 
machine plots. 


Input Data Printout 

The card images of all the input data - configuration description and plot cards - 
are printed. Cards for a sample case input deck are listed in table I. 

Cross-Section Plot Printout 

The points of intersection of the cutting plane and the given body may be printed as 
specified on the cross-section plot card. The actual points in the original coordinate 
system, the actual points in the scaled coordinate system, the rotated and projected points 
in the original coordinate system, and the rotated and projected points in the scaled 
coordinate system may be chosen for printing or printing of the intersection points may 
be omitted entirely. Since the plane of intersection is always transformed into the 
YZ -plane, the rotated x-coordinates should be all zeros. The equation of the plane is 
also printed. 


Plot Vector File 

A plot vector file is generated during execution of both plot options - three- 
dimensional orthographic or cross sectional. The plot vector file can be postprocessed 
on the same computer run or at a later time for the desired plotting device. 

MACHINE SETUP 

This program was written in FORTRAN Version IV for Control Data series 6000 
computer systems with the Scope 3 operating system and library tape as modified for the 
Langley computer facility. Tape unit 5 is used for input, unit 6 for output, and units 7 
and 10 for intermediate storage. Approximately 60000 octal locations of core storage 
are required and the processing of information for one plot is less than 1 minute of com- 
puter time. 


OPERATIONAL DETAILS 

The graphic output system at Langley Research Center is in two parts: (1) a device 
independent graphic language which produces a plot vector file containing only plotting 
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commands, and (2) a set of postprocessors which format the output of the graphic language 
to a particular graphic device or devices. 

Subroutines PSEUDO, NFRAME, NOTATE, and LINE are the basic subroutines used 
from the graphic software package. Subroutine PSEUDO causes the necessary parameters 
and linkage to be set up to output a device independent plot vector file for postprocessing 
for a particular plotting device. Subroutine NFRAME provides a means of executing a 
frame advance movement. Alphanumeric information for annotation and labeling is drawn 
by subroutine NOTATE. Subroutine LINE draws a continuous line through a set of succes- 
sive data points where the minimum values and scale factors are stored at the end of the 
data arrays. 


CONCLUDING REMARKS 


A computer program has been written to define mathematically an arbitrary curved 
surface in surface -patch -equation form. Although the program is oriented toward air- 
craft configurations, it can be used to model mathematically any three-dimensional object 
by using an alternate data input format. Parametric cubic spline curves are fitted to the 
input data points to define the boundary -curve slopes. These slopes and the corner points 
are the necessary components of the surface patch equations. Program options include 
the application of three-dimensional rotation equations directly to the patch equations for 
plotting the surface at any desired viewing angle and the solution of a number of patches 
and a plane for generating cross-section or contour plots of an object or surface. Output 
from this program has been used to drive Calcomp, Gerber, and Varian plotters. The 
program has also been used for on-line display on a cathode -ray-tube device. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., February 28, 1975. 
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PARAMETRIC CUBIC SPLINE SPACE CURVES 


The derivatives in u and w at the corners of the surface patches are required 
for the surface equation. The grid points which determine the boundary curves are fitted 
in the w-direction and then in the u-direction for computation of the derivatives. A 
parametric cubic spline curve fit technique by Timothy E. Johnson of Massachusetts 
Institute of Technology is used. 

Parametric cubics are the lowest order polynomial with the property of being able 
to twist through space. The spline curve exhibits the usual meaning of smoothness by 
minimizing curvature. Parametric curves are not sensitive to infinite slopes. 

For purposes of completeness a brief description of the method follows as well as 
an explanation of the calling sequence for Subroutine SPFIT which utilizes this method. 

A series of adjacent polynomial segments between each pair of given points is 
used to represent the curve. 

The component cubic parametric equations are 
X(t) = A x t 3 + B x t 2 + C x t + D x 

Y(t) = A y t 3 + B y t 2 + C y t + Dy 

Z(t) = A z t 3 + B z t 2 + C z t + D z 

where all operations are performed once on each component equation. 

The component equations will be represented by 

P(t) = At 3 + Bt 2 + Ct + D (1) 

The tangent parametric slopes at each point define the coefficients of the cubic seg- 
ments and must be determined so that the tangent slopes give a smooth curve. 

The given points Pp P 2 , . . ., P n have the corresponding tangent slopes 

p'. p* . . p ? and the independent variable t varies so that 0 s t S L where L 
1 ' z n. 

is the chord length between the given points. (See sketch (f).) 
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The boundary conditions for one cubic segment P(t) can be written as 
P(0) = Pj 

P ( L l) = P 2 


dP(t) _ p* 
hT“7 - P 1 


dP(t) 


dt 


t=Li 



The coefficients of equation (1) can be written in terms of the end points and para- 
metric slopes for a given segment 


P(0) = D 1 = Pj 


_ c = p* 

dt 1 1 


t=0 


and 


P(L X ) = AL 1 3 + BLj 2 + PjLj + Pj = P 2 

■ dP ^ = 3AL, 2 + 2BL-, + Pi, = P' 
dt 1 112 

at t=L 1 

Solving the last two equations yields 
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A _ 2 ( P 1 - P 2) , 
A i — n 


Vfi. 
l i 2 w 2 


B 





The coefficients for each adjoining segment are found in the same manner. 

The adjacent cubic segments can be related by setting the second derivatives equal 
at the common points. 

For equal second derivatives at P 2 


d 2 P(t) _ d 2 P(t) 


dt: 


t=Li 


dt: 


t=0 


6A x Lj + 2Bj = 2B 2 (2) 

Substituting the expressions for A and B into equation (2) and collecting terms give 
in terms of the unknown slopes: 


L 2 P 1 + 2 ( L 2 + L l) P 2 + L 1 P 3 " ( P 3 " P 2) 



This equation may be written over all segments in matrix notation: 

L 2 2 ( L 1 + L 2) L 1 

O 

— 


Pi 


3 / L l L 2[ L l 2 ( P 3 ' P 2) + L 2 2 ( P 2 * P l)) 

L 3 2 ( L 2 + L 3) 

L 2 



P 2 


3 / L 2 L 3 [ L 2 2 ( P 4- p 3) +L 3 2 ( P 3- P 2)] 

L 4 

2(L 3+ L 4 ) 

L 3 

X 

P ’ 3 



o 






/ c 2 2 / n 





P n 


3 / L n-2 L n-l[ L n-2( P n * P n-l) + L „-l( P n-l ' P n- 2 )J 
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Since there are two more unknowns than there are equations, the second derivatives at 
Pj and P^ are made equal to zero. The Thomas Algorithm, which is equivalent to 
Gaussian Elimination without pivoting, is used for solving the tridiagonal matrix. The 
cubic coefficients are normalized so that the parametric space for each curve segment 
varies from 0 to 1: 

A’ = AL 3 

B' = BL 2 

C’ = CL 



APPENDIX B 


BICUBIC SURFACE PATCHES 


This appendix is a brief description of the surface patch method used in pro- 
gram D3400. 

The x-, y-, and z-coordinates of surface points are functions of two variables 
u and w: 


X = f(u,w) 

Y = g(u,w) 

Z = h(u,w) 

V = (X Y Z) 

The vector V(u,w) is a function of the two variables, u and w, which take on only the 
values between 0.0 and 1.0. 

The method builds a surface by joining surface "patches." A surface patch can be 
thought of as a portion of a surface bounded by four space curves: (0,w), (l,w), (u,0), 

and (u,l). 

The surface patch equation used is 
V(u,w) = UMBMM 

where 

U=(u 3 u 2 u l| 

W = (w 3 w^ w ^ 


V(0,0) 

V(0,1) 

v w (0,0) 

V W (0,1) 

v(l,0) 

V(l,l) 

V W (1,0) 

V w (l,l) 

V u (0,0) 

v u (0,l) 

VuwW.O) 

Vuw* 0 ’ 1 ) 

V U (1,0) 

V u (l,l) 

V UW (1,0) 

v w <m> 
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This 4x4 matrix is called a boundary matrix as it contains only geometric properties 
from the boundary curves: corner coordinates, corner slopes, and corner twists. (The 
corner twists are all made equal to zero in the present application.) These quantities 
are constants and grouped systematically in the matrix. 


The blending function matrix M provides a blending effect from one surface patch 
into an adjoining patch 



Since the matrix product MBM t is constant, the pre- and post-multiplications are 
performed and the equation stored as S = MBM* for further computer processing. It is 
a simple matter, then, for writing an expression for u with w held fixed or u held 
fixed and w varying. 

For a detailed development of the surface patch equation, see reference 1. Refer- 
ence 2 presents the bicubic form of the surface patch equation as used in this report. 
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DESCRIPTION OF FILE AND METHOD OF STORAGE 
FOR SURFACE PATCH EQUATIONS 


Program D3400 computes and stores by columns of the surface grid the patch 
matrix equations using the previously stored lines of grid points. A disk file is used for 
temporary storage of the elements of all the patch matrix equations. 

The elements and certain other information are written in binary on disk unit 7 in 
the following format. FORTRAN names from the given program are used for illustration. 


Record 


Name 
ABC(8) 

XMIN Minimum X of the entire body 

XMAX Maximum X of the entire body 

YMIN Minimum Y of the entire body 

YMAX Maximum Y of the entire body 

ZMIN Minimum Z of the entire body 

ZMAX Maximum Z of the entire body 

NOBJ Total number of objects (or components) to form the 

complete body 

(for NOBJ 1 ) 

NSURF Number of surfaces in object 

Ml Not used 

M2(2) , Name of object (for printing only) 

M3 Not used 

M4 Not used 

(for NSURFj) 

Ni Number of rowwise patches (NCOL - 1) 

Number of columnwise patches (NROW - 1) 

N3 Not used 

N4 Not used 

N5 Not used 

PATCH (4,4,3) 48 elements of a patch equation are written by columns 

as one record 

(Record 5 is repeated for N1 x N2 patches taken by columns over the 
surface grid) 

(Records 4, 5, and so forth are repeated for each NSURF from record 3) 

(Repeat record 3 for each NOBJ as given in record 2 and repeat records 4, 
5, and so forth as required) 


(for NOBJ 1 ) 
NSURF 
Ml 
M2(2) 

M3 


(for NSURFj) 
Nl 


PATCH (4,4,3) 


Purpose 


Identification 
Minimum X 


Maximum 

Minimum 

Maximum 

Minimum 

Maximum 
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ORTHOGRAPHIC PROJECTIONS USING 
SURFACE PATCH EQUATIONS 

The orthographic projections illustrated in this report are created by rotating each 
patch element to the desired viewing angle, transforming the rotated patch into a coordi- 
nate system in the plane of the paper, and computing actual coordinates from the patch 
equation. The body coordinate system is coincident with the fixed system in the plane 
of the paper when all the rotation angles are zero; for example, the configuration X-axis 
and Y-axis would coincide with the paper for plots in the X’Y' paper plane. 

The rotations of the body and its coordinate system to give a desired viewing angle 
are specified by angles of yaw, pitch, and roll (xj/, 9, and <p), shown in sketch (g). 


if 



Sketch (g) 


The equations used to transform the patch equations defining the body with a set of 
rotation equations (xj/, 6, and </>) into the desired paper plane are 

Sx = S x (cos 9 cos xj/) + S y (-sin xj/ cos </> + sin 9 cos x p sin <p) 

+ S z (sin xj/ sin <fi + sin 9 cos xj/ cos (p) 

Sy = S x (cos 9 sin xj/) + Sy(cos xp cos (f> + sin 9 sin xj/ sin </>) 

+ S z (-cos xj/ sin <p + sin 9 sin xj/ cos 0) 
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s z = S x (- sin o) + S y (cos e sin 0) + S z (cos 9 cos 0 ) 


Each element of the patch equation is rotated to the desired viewing angle and then 
transformed into a coordinate system in the plane of the paper. Only two of the preceding 
equations, determined by the desired paper plane, are used and result in a two-component 
patch equation. 

The body surface plots may be enriched to any desired degree. By assigning u a 
value between 0.0 and 1.0, we have two cubic equations in w which may be easily solved 
by varying w from 0.0 to 1.0 for the rotated and projected coordinate values of a line 
across the patch. In a similar manner, w may be assigned a value from 0.0 to 1.0 
resulting in cubic equations in u. By varying u from 0.0 to 1.0, coordinate values are 
generated across the patch in the u-direction. 

An optional hidden-line test is incorporated into the program which may be used to 
provide the capability of deleting most elements on the surface of the configuration which 
would not be seen by a viewer. No provision is made for deleting portions of an element 
or components hidden by other components. 


The surface vector normal to the paper plane is computed at each interpolated point. 
(See ref. 1.) If the surface vector is positive, the computed point faces the viewer and is 


visible; 

if the surface vector is negative 

the computed point is not visible and is not 

plotted. 

With 






£ 

II 

Y u 


9Y 

3u 

azl 

auj 


J* 

Cj 

ii 

£ 

Y w 

z w) =(f 

9Y 

aw 

az] 

awj 


For one surface, the surface normal vector is 


N=UXW = [J X J y 

where 
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z u 

Xu 

z w 

Xw 

Xu 

Y u 

Xw 

Y w 


The visibility of the point under consideration is determined by evaluating only the 
Jacobian outwardly normal to the paper plane at the point. 
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CROSS-SECTION OR CONTOUR PLOTS USING 
SURFACE PATCH EQUATIONS 


By using both the surface patch equation in the form 
V(u,w) = USW T 

where the components of the 4 x 4 matrix S are 0S X S y S0 (see appendix B) 
and the equation of a plane 

ax + by + cz - d = 0.0 

and equation of the intersection of the plane and the patch surface may be written. Sub- 
stitute x = US X W T , y = US y W T , and z = US Z W T in the preceding equation so that 

U(aS x + bSy + cSjW T - d = 0.0 

The matrix 00 = QiS x + bSy + cS0 is composed of constant elements and may be eval- 
uated for an equation in the two variables u and w 

UGW T = d 

By assigning values to w, a series of cubic equations in u are solved for points on the 
intersection curve of the surface with the plane or values may be assigned to u for 
cubic equations in w to solve for the intersection curve. (See ref. 1.) Only one solution 
for a curve of intersection is allowed. 

The points of the curve of intersection are then rotated and translated so that the 
plane of intersection coincides with the YZ-plane of the paper. The method for rotation 
and translation follows. 

Sketch (h) illustrates the transformation system established by the 3 input points 
used to define the cutting plane. 
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z 



The translated and rotated coordinate system used for plotting is shown in sketch (i). 


Sketch (i) 

Sketch (j) illustrates the choosing of three points on a cutting plane when the method 
of defining the plane is by the input of an X-intercept, a roll angle 6 , and a Mach number. 

XGF One-quarter of Mach cone 

EFGHIJ Lies in plane parallel to YZ-plane 
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JXH Mach plane, tangent to cone along IX 

IE Projection onto plane JHE of normal to Mach plane at I 

To locate three points on the plane: 


Point 

X 

Y 

Z 

P 1 

X 

0.0 

0.0 

P 2 

X - c 

0.0 

1.0 

Ps 

X - B 

1.0 

0.0 


where the equation of the plane is expressed as AX + BY + CZ = D with 
/3 = M^ - 1, A = 1.0, B = -j3 cos 9, C = -j3 sin 9, and D = X 


z 



By using the three input coordinates that define the cutting plane, two diagonal vec- 
tors may be formed with the components: 
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T 1X " x 3 " X 1 
T 2X * x 3 - *2 


T lY =y 3->'l 
T 2Y “ s 3 - y 2 



The cross products of the diagonal vectors give 


N X " T 2Y T 1Z " T 1Y T 2Z 


N Y “ T 1X T 2Z " T 2X T 1Z 


N Z ~ T 2X T 1Y “ T 1X T 2Y 


I 2 o 2 

which when divided by N = \l N x + Ny + N z yield the unit normal vector n of the 
cutting plane 


n NX 
x ¥ 



n.1-7 “ 



To determine the coordinate system for plotting, two unit vectors lying in the cutting 
plane are needed 


v 





where 

T = \/ T ix 2 + t iy 2 + T 1Z 2 

and 
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V 2X = n Y V lZ ‘ n Z V lY 
V 2Y = n Z V lX " n X V lZ 
V 2Z = n X V lY " n Y V lX 


The transformation matrix becomes 



11 

a 12 

a 13 

21 

a 22 

^3 

31 

a 32 

a 33 


where 


a H = n x 

a 12 = n Y 

a 13 = n Z 

a 21 " V 1X 

a 22 = V 1Y 

a 23 = V 1Z 

a 31 = V 2X 

a 32 = V 2Y 

a 33 = V 2Z 


The origin of the reference coordinate system is transferred to the first point given 
in defining the cutting plane 


x 0 


' x i 

yo 

-ra 

y i 

A 


_ z i 


The transformation from the reference to the paper coordinate system is accom- 
plished by: 
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TABLE I.- SAMPLE CASE INPUT CARDS 

SIMPLE CAMBERED CIRCULAR BODY 


1 1 

-1 1 

1 

6 5 

1 5 

6 

9494. 

0 . 

20 . 

50. 

70. 

100 . 


82 • 30 

5.05 

0 . 

180. ICO 


I 14.1999.90 

-.45 

142.351 


157.98 

19.80 

-1 .85 

98.570 



239.18 

46.20 

-2.80 

36.719 



269.23 

59.40 

-4.30 

15.670 



282.00 

66*00 

-4.40 

7.400 



3.6 

2.75 

-3*45 

- 6.8 

-9.4 


0 . 

1 *35 

-1 *2 

-3.45 

- 6.8 


0 . 

• 72 

• 0875 

-.7825 

-2.173 


0 . 

• 248 

• 31 1 

• 2995 

• 2385 


0 . 

-.043 

-.090 

-.110 

-• 1224 


0 . 

-.0325 

-.075 

-. 1 00 

-• 1 324 


0 . 

1 .069 

1 .518 

1.451 

0 . 


0 . 

• 889 

1 .272 

1.136 

0 . 


0 . 

• 086 

1 .294 

1 .087 

0 . 


0 . 

• 880 

1 .375 

1.155 

0 . 


0 . 

• 880 

1 .375 

1.155 

0 . 


0 . 

• 880 

1 .375 

1.155 

0 . 


0 . 

70. 

130. 

200 . 

260. 

312. 

7.4 

7.4 

1 .25 

-7.45 

- 10.2 

- 10.2 

0 . 

96. 

98. 

79. 

62. 

0 . 

241 .0 

31.75 

-3.60 




0 . 

12 . 

24. 

34.5 



2.292 

2.79 

3.08 

3. 1 



252.0 

47.0 

-2.95 

35.3 

285.36 

47.0 

0 . 

30. 

50. 

70. 

100 . 


0 . 

• 76 

• 977 

• 927 

0 . 



1 


X z 
x z 
X z 

2 

1 0. 
-45. 
-45. 
OUT-45. 

4. 

10 . 
1 o. 
10 . 

8 . 

-30. 

-30. 

-30. 




200. 

0. 

o. 

200 . 

0 . 

1 • 

200 . 

200. 

0. 

0. 

200 . 

0 . 

1 • 

200 . 

-10. 

2.5 

0. 

0 . 

0 . 

45. 

1 .2 


-10. 

5. 

0. 

0 . 

55. 

45. 

1 .2 



14 15 

REFA 

XAF 5 

WAFORG 1 
WAFORG 3 
WAFORG 5 
WAFORG 9 
WAFORG 1 1 
WAFORG 1 2 
TZORD 1 
TZORD 3 
TZORD 5 
TZORD 9 
TZORD 1 1 
TZORD 12 
WAFORD 1 
WAFORD 3 
WAFORD 5 
WAFORD 9 
WAFORD1 1 
WAFORDl 2 
XFUS 6 
ZFUS 6 
AFUS 6 
PODORG 2 
XPOD 
POOR 

FINORG 2 

XF IN 

F I NORD 



15. 

VU3 


2 




15. 

ORT 


2 




15. 

ORT 5 

6 

2 




15. 

ORT 5 

6 

2 

1 


1 • 

0 . 

.30 

0 . 

16.PNT 



1 • 

0 . 

.30 

20 . 

16. PNT 

8 

8 



• 1 

0 . 

1 6 • A NG 

8 

8202 



• 1 

20 . 

1 6 • A NG 

8 

8482 
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Figure 3.- Alternate input from patch definition. 
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Figure 6.- Optional input cross sections with constant reference. 
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